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Chapter  I  INTRODUCTION  AND  LITERATURE  SURVEY 
A.  Introduction 

An  important  development  in  recant  years  has  been  the 
evolutions  of  chemical  and  physical  techniques  to  produce 
highly  ordered  polymers.  Some  of  the  resulting  materials 
exceed  all  others  in  their  ratio  of  tensile  strength  to 
density  ( s / p ) .  Because  of  their  potential  for  applications 
capitalizing  on  their  large  s/a  value  the  major  portion  of 
research  on  ordered  polymers  has  dealt  with  mechanical 
properties  and  with  methods  of  synthesis  or  production. 
Other  subjects,  such  as  transport  properties ,  have  received 
relatively  less  attention.  As  pointed  out  by  Barker,  Tsai, 
and  Wiliency,^  transport  properties  such  as  electrical  con¬ 
ductivity  and  the  diffusion  of  various  molecular  species 
in  polymers  not  only  are  of  interest  in  themselves  but  also 
provide  useful  probes  to  gain  a  better  understanding  of  the 
polymer's  properties  in  general.  Some  of  the  most  recently 
discovered  types  of  ordered  polymers  are  the  extended  chain 
(rigid  rod)  polymers.  The  main  material  studied  in  this 
thesis  is  the  trans-isomer  of  poly  (paraphenvlene- 
benzobisthiazole)  (PPBT) . 


,  r  -  -  - 


h.,w. 


The  available  PPBT  samoles  were  of  very  small  dimensions 


(the  fibers  were  only  10  um  -  4  0  ym  in  diameter  and  the 
films,  only  10  pm  -  30  ym  thick  and  5  mm  wide).  Furthermore, 
the  available  samples  exhibited  considerable  internal  and 
external  imperfections  which  may  lead  to  certain  difficul¬ 
ties.  At  this  point  it  is  useful  to  summarize  some  of  the 
problem  areas: 

(i)  Difficult  manipulation  due  to  the  small  sizes 
of  samples. 

(ii)  Variable  diameter  along  a  fiber  and  uneven 
surface  and  thickness  of  a  film. 

(iii)  Small  total  quantity  of  compounds  sorbed  by 
the  polymer. 

Regarding  difficulty  (i)  the  experiments  are  hard  to 
perform  not  only  because  of  the  small  size  but  also  due  to 
a  tendency  of  the  samples  to  fibrilate.  Difficulty  (ii) 
will  induce  experimental  error,  because  diffusion  depends 
strongly  on  the  thickness.  Difficulty  (iii)  requires  that 
the  experiments  be  carried  out  with  highly  sensitive 
instruments . 

Since  it  so  ubiquitous  and  plays  such  an  important 
role  with  regard  to  the  electrical  properties  of  polymers, 
one  of  the  major  diffusants  studied  was  water.  It  exists 

and  affects 


in  the  normal  atmosohere 


the  orcoerties 


of  PPBT. 


3 


The  reaction  between  PP3T  and  Ho0  was  assumed  to  be  a 
2 

prctonation  effect.  From  an  observed  shirt  of  the  C=N 
stretching  band  in  laser  Raman  spectra  from  1481  cm  ~ 
(solid)  to  1515  cm  ~  (in  methane  sulfonic  acid  (MSA) )  the 
proconation  process  in  MSA  was  deduced  to  be  expressed  by 
the  following  equation  where  ?  denotes  the  neutral  PPBT 
chain : 

P  +  nCH-HSO,  t  H  +  nCH-.SO~. 

3  3  n  3  3 

Then  for  water,  the  analogous  reaction  is: 

P  +  nH„0  t  H  ?n+  +  nOH~. 
z  n 

Because  of  this,  moisture  affects  the  electrical  conductivi 
of  PPBT  ir.  a  strong  way.  This  has  been  extensively  studied 
by  D.  Y.  Chen.3 

In  order  to  study  the  diffusion  of  a  given  species 
through  PP3T ,  three  methods  have  been  developed:  1.  X-ray 
energy  dispersive  analysis.  2.  Weighing  method  with  a 
Cahn  electrobalance,  and  3.  Sorption  and  Desorption  for  a 
sample  suspended  on  a  quartz  spring  in  an  evacuable  chamber 
Besides  water,  ethanol,  and  benzene  also  have  been 
studied.  Benzene  v;as  chosen  because  its  size,  shape,  and 
electronic  structure  relative  to  those  properties  of  the 
mer  make  it  of  special  interest.  The  polymer  has  been 


studied  under  different  temperatures  so  that  the  activation 
energy  E^,  solubility  s  and  the  preexponential  factor 
have  been  obtained.  The  data  from  experiments  in  the  Cahn 
electrobalance  were  not  consistent  with  those  from  experi¬ 
ments  using  the  quarts  spring  balance.  Efforts  to  resolve 
this  conflict  led  to  the  discovery  of  a  useful  model  based 
on  a  distribution  of  microcavities.  A  molecular  medal  in¬ 
volving  the  partial  rotation  of  chain  segments  has  been 
developed  to  account  for  the  activation  energy.  This  model 
is  consistent  with  the  chain  packing  configuration  as  in¬ 
ferred  by  X-ray  diffraction. 

When  PPBT  was  doped  using  ionic  compounds  (nitrates 
and  chlcrid.es)  to  provide  trapping  sites,  the  diffusion 
coefficients  of  Ho0  in  doped  PPBT  vs  the  charge  to  radius 
ratio,  q/r  of  the  cations  was  found  to  be  linear. 

Finally,  the  parallel  diffusion  in  PPBT  has  been  con¬ 
sidered  from  the  view  points  of  fibrilar  morphology  and 
crystal  structure. 

For  a  thorough  understanding  of  the  diffusion  processes 
in  PP5T,  we  must  be  able  to  represent  the  mass  flow  by 
means  of  some  equations  based  on  different  experimental 
methods  and  theoretical  concepts.  The  following  pages  of 
background  discussion  will  serve  this  end. 
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3.  Literature  Survey 

1.  Pick's  Laws.  Approximately  a  century  ago  (in  1889) 

Pick  developed  a  continuum  approach  to  the  phenomenon  of 
diffusion  which  still  provides  the  framework  for  discussion 
of  the  principles. 

Assume  diffusion  only  in  the  +x  direction,  then  consider 
a  thin  layer  of  diffusant  between  x  and  x+dx.  The  number 
of  molecules  of  diffusant  (of  type  j)  per  unit  volume  at  any 
time  t  is  a  function  only  of  x.  Let  Cj(x,t)  denote  this 
concentration.  The  diffusive  flux  Jh  of  molecules  across  a 
perpendicular  plane  at  x  is  the  net  number  of  molecules  of 
j  passing  through  unit  area  of  the  plane  in  the  positive  x 
direction  in  unit  time. 

At  the  molecular  level  diffusion  is  a  statistical 

phenomenon  which  can  be  related  in  a  useful  way  to  changes 

in  the  chemical  potential  of  the  diffusant  in  the 

diffusate.  The  chemical  potential  may  be  viewed  as  the 

statistical  analogy  of  the  potential  energy  function  in 

classical  mechanics.  Thus  the  gradient  of  Uj  in  a  given 

direction  x  corresponds  to  the  averaae  force  on  a  diffusant 
a 

molecule  due  to  interactions  with  its  molecular  surroundings. 

,1  a  .  . 

f  A  =  -  •  rr-  ,  ( a  =  x  ,  y ,  z )  ( 1 ) 


oar 


is  included  to 


The  factor  of  l/Avogadrc's  nurru 
convert  the  usual  units  for  chemical  potential  (J/ir.ol)  into 
energy  per  molecule,  so  that  f_.  is  the  i-coaponant  of  for. 
per  molecule  of  type  j.  This  average  force  causes  a 
statistical  drift  of  tne  drffusarc  with  a  snood  u.  given 

1  ,  a 


u  . 
J 


(2) 


where  w.  is  the  molecular  mobility  of  species  j.  The  net 
diffusive  flux  densitv  is 


"j  .a  nj'i«u 

with  C.  =  n./N„  relating  the 
J  J  A  y 

that  n.  in  molecules /m^ .  The 
1 

to  the  concentration  through 
relation 


-C  o  -  (  3  u  . 
3  3  3 


/Ox, ) N, 

‘.C  r-» 


-1 


(3) 


3 

concentration  C.  in  mol/rr,  to 
chemical  potential  is  related 
the  familiar  thermodynamical 


u .  =  u°  +  RTln (y -C . )  (4) 

3  3  3  3 

from  which  the  usual  relation. 


-D  (  3C  .  / Ox  ) 

3  /  J  i 


(5) 


Known  as  Fick's  first  law  is  seen  to  follow  as  a  restricted 

soecial  case  (when  the  activity  coefficient  v.  is  not 

,  -1 
-i 

dependent  on  and  when  there  is  no  thermal  gradient  and 
J 

no  interphase  boundary) . 


.  S  V 
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J  (x)  =  (omitting  the  subscript  j) 

6  X 


In  the  absence  of  chemical  reactions,  etc.  the  conserve t 
of  component  j  may  be  expressed:  Accumulation  rate  = 
influx  rate  -  outflow  rate 


(Adx,3C  _  A [  J  j x )  -  J(x  +  dx)ldx 
'  X;3t  ax 


_  _  _  a  ,  „  c;  l.  . 

3x  3x  lu~3x' 

=  Co— %  +  ^(Cu) 

3x2  "X  "X 


or  for  the  conditions  appropriate  cor  Pick’s  first  law,  we 
have  Pick’s  second  law 

3C.  3  2C . 

-K-r  =  (o.kT)  - ^  ,  D.  =  v).kT  (B) 

dt  3  w-  0  J 


2.  Experimental  Methods. 

a.  X-Ray  Energy  dispersive  analysis.  This  me  shod  'was 
developed  by  Barker,  Reynolds  and  Malhotra  ,  and  was  successful! 
applied  to  llg  diffusion  in  Ag^S.n  Amalgams. 

The  solution  to  the  on-dimensional  diffusion  equation 
for  a  semi- infinite  solid  is  (Appendix  A) 


=  C  orfc(~ ] 
°  >'4Dt 


erf  (— — ) 
V  4  D  t 


3 


prior  to  the  establishment  or  a  steady  state,  both  the  rare 
of  flow  and  the  concentration  at  any  point  within  the  film 
vary  with  time.  If  the  diffusion  coefficient  is  constant, 
the  film  is  initially  completely  free  of  diffusant,  and 
diffusant  is  continually  removed  from  she  low  concentration 
side,  the  amount  of  diffusant  Cv  which  passes  through  the 
film  of  thickness  b  in  time  t  is  given  by 


t  =  +  1  _  4,  y  1-y  '  e  ~  '  b 


nw  2 

-i  i71  ~D t- 


AC~b  , 
2  b 


2  L  .2 

a  n=i  n 


(li) 


As  t  =°,  the  steady  state  is  approached  and  the 
exponential  terms  become  negligibly  small,  so  that  the 
graph  of  Qt  against  t  tends  to  the  line 


Qt/A  =  (DC2/b)  (t  -  b“/6D) 


(11; 


which  has  an  intercept  x,  on  the  t-axis,  given  by 


x  =  b  /6D. 


(12) 


"heref  ore 


D  =  b  2  /  6  - 


(13) 


The  mathematical  details  are  shown  in  Appendix  B. 


c.  Sorotion  and  Dcsorotion  Method.  If  the  concen¬ 


tration  just  within  the  surface  of  a  film  exposed  to 


10 


diffusant  is  maintained  constant,  then  the  amount  of 
diffusant  ,  taken  up  by  the  film,  of  thickness  b  in  a  ti.r.. 
t,  is  given  by  the  equation 


o  1 

m.  /m  =  i  -  4  y  — 

tt  n=G  (2n+l) 


r  D(2n+i)2m2t1 
exp  [ - - - - - ] 


Aooendix  C 


:s  the  details. 


M  is  the  equilibrium  sorotion  attained  thsorsticallv 
after  infinite  time.  Assume  t,  is  the  time  at  M.  /M  = 
the  above  equation  can  be  written  as 


Vb  = 


Thus ,  we  have 


1  m  1  q 

-y-  Inljr-  -  1 

TT 2  D  10  9 


D  =  0.04919  b  / 1 j 


The  time  t,  can  be  found  experimentally  and  thus  D  can  be 
i 

calculated.  Another  form  of  equation  (14)  is 


2  k  h  7.  n 

Mfc/M  =  4  (Dt/b  i‘(l/r  +  2  l  (-1)  nierfc  (nb/2  (Dt)  ')  )  (17) 

u  n=0 


2  2 

ierf  (nb/2  (Dt)  '2)  =  (l/,7)e_n  b  /4Dt 


-  (nb/2  (Dt)  Merfc  (nb/2  (Dt)  2) 
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•Ths  value  of  D  can  be  deduced  from  an  observation  of  the 
initial  gradient  of  a  graph  of  M^/M^  as  a  function  of 

7  k 

(t/b“)  2. 

This  kind  of  analysis  has  been  used  by  Garrett  and 

7 

Park  (1967)  (Garrett,  I960)'  for  kinetic  and  equilibrium 

3 

sorption  measurements.  Others,  for  example  Baughan,  and 

9 

Carpenter  and  Twiss  were  also  working  on  these  kind  of 
methods . 

The  equations  given  above  are  for  films.  The  same 
approach  has  been  applied  to  fibers  and  spherical  particles. 
For  fibers  the  diffusion  coefficient  D  was  determined  from 
the  slope  of  the  linear  relationship  obtained  when  plotting 

j. 

(Mt/Mca)  2  against  t  using  the  following  equation 

D  =  Ur2/16t)  (M./M  )  ^  (19) 

t 

d.  An  Optical  Method.^  This  method  was  applied  to 

diffusion  in  which  the  diffusant  made  a  visible  border  in 

polymers.  There  are  two  known  situations  in  which  sharp 

borders  occur.  The  first  is  the  most  widely  known.  Namely 

swelling  of  the  polymer  by  the  diffusant  creates  an  advancing 

zone  between  the  swollen  phase  and  the  unswollen  phase.  To 

gee  such  an  effect  D  ..  >  D  The  second 

swollen  unswollen 

situation  occurs  when  there  are  colored  trapping  sites  which 


change  color  or  refractive  index  when  a  diffusant  molecule 


is  trapped.  J.  J .  Hermans  proposed,  a  mechanism  m  which 

the  boundary  line  indicates  the  abrupt  change  from  "all 

holes  filled"  to  all  holes  empty.  Figure  (1)  gives 

diagramatic  representation  to  the  situation. 

i  o 

Barker  used  this  idea  and  derived  an  equation  that  is 
very  practical.  He  assumed  the  shaded  area  in  Figure  l.b. 
has  approximately  the  area  of  a  triangle,  then  from  equation 


J  = 


He  obtained  the  equation 


-DdC/dx  =  n  (d,',/dt)  +  (C  / 2)  (d^/dt)  a  D Cj? 


Mere,  J  is  the  amount  of  diffusant  entering  the  polymer  in 

unit  time.  corresponds  to  the  increased  concentration  of 

free  molecules,  J ^  represents  the  flux  going  into  traps 

between  £  and  £+d£,  and  n^  is  the  concentration  of  traps. 

The  concentration  maintained  just  inside  the  surface  of  the 

ooivmer  is  denoted  bv  C  . 

s 

If  C  and  n~  are  constants 


D  ^  U“/2t)  (nQ/Cs  +  h) 


later  on,  Barker  found  from  experiments  for  C>2  diffusion  in 
polycaroonate  that  better  agreement  is  obtained  by  replacing 
the  h  by  0.32. 


e.  Other  Methods.  Certain  other  methods  can  detect 

concentrations  at  a  particular  distance  from  a  boundary, 

for  example  certain  refractive  index  techniques,  radiation 

1  2 

a o sorption  metnods  ( spec trophctome trie  methods),  etc. 
According  to  the  equation 

C  (x  ,  t )  =  (C.,/ (aDt)i5)exp<-x2/4Dt)  (22) 

O 

which  represents  the  concentre  tier,  profile  at  time  t 
subsequent  to  an  initial  highly  localized  ccn.cen tra  cion  on 
a  planar  interface.  On  plotting  lcgC(x,t)  vs.  xz,  a 
straight  line  having  the  slope  1/  (4Dt)  should  be  obtained 
from  which  D  may  be  calculated  if  t  is  brown . 

2.  Seme  Other  effects  Related  to  Diffusion. 

a.  Relations  between  D,  s,  and  ?.  The  diffusive 
transport  through  homogeneous,  noncorous  polymer  membranes 
was  qui  .e  generally  treated  as  if  the  membrane  were  ideal 
so  that  the  flow  could  be  described  by  Pick's  laws  based  on 
the  concentration  of  the  permeant.  In  such  a  case  the 
current  density  J  through  the  membrane  is  the  product  of 
the  negative  concentration  gradient  and  the  diffusion  co¬ 
efficient  (Pick's  first  law  of  diffusion) 


J 


-ODC/  5.x 
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In  a  planar  mambrane  of  uniform  material  properries, 
the  steady  state  concentration  gradient:  is  constant  so  that 
the  flux  per  unit  area  of  the  membrane  is  given  by 

J  =  -DbC/ix  =  D(C1  -  C2)/b.  (22) 

where  the  membrane  has  a  thickness  b  which  is  uniform  and 
D  is  independent  of  C.  The  subscript  1  relates  to  the 
upstream  and  the  subscript  2  to  the  downstream  boundary  of 
the  membrane. 

In  case  of  an  "ideal  gas-membrane  pair"  one  has  the 
simple  relationship  for  an  isothermal  sorption  process 

C  =  sp  (24) 

which  defines  s  and  which  is  just  Henry's  law  if  s,  the 
sorption  coefficient,  is  a  constant  independent  of  pressure 
p,  then  from  equation  (23), 

J  =  -Dsip/Ax  =  Ds(p^  -  p2 )  /b  =  -  ?,\p/Ax  (25) 

so  that  the  constant  permeability  coefficient  is  given  by 

P  =  Ds ,  (26) 

although  there  u.ju.i  *_Iy  is  an  extra  constant  factor  to 
handle  various  mixtures  of  unit  systems. 


Dependent 


on  Co 


Usually  organic  vapor  diffusion  is  concentration-dependent . 
Crank  introduced  a  mean  value  5,  which  is  the  average  of  the 
diffusion  coefficients  for  a  number  of  different:  surface  con¬ 
centrations.  Calculations  have  shown  that  for  any  one 
experiment,  5  provides  a  reasonable  approximation  to 

Cf 

1/C*  J  D  dC  (27 

"  0 

where  0  to  Cc  is  the  concentration  range  existing  in  the 


sheet  during  that  experiment, 


If  D  is  calculated  from  th 


initial  gradient  for  sorption  and  from  the  gradient  for 
desorption,  then  (Dg  +  D^)  is  a  better  approximation  to 


(1/Cf)  j  D  dC 


c.  Temperature  Effect.  Under  ideal  conditions, 
diffusion  follows  the  Arrhenius  eauation 


-Ed/kT 


D  =  D0a 


over  sufficiently  restricted  temperature  ranges.  The  pre- 
exoonential  factor  D„  includes  entropic  terms. 

a.  Gas  or  Vapor  Properties  that  Affect  Diffusion . 
Usually,  larger  gas  molecules  need  1  iraer  holes  in  the  rolyrv 
for  diffusion  to  proceed.  Accordingly,  the  smaller  the 
molecules  are,  the  smaller  the  activation  energy  will  be  and 


con.32quer.tIy  the  larger  the  diffusivity  will  be.  If  the 

molecular  size  is  large,  the  effects  of  molecular  shape  and 

flexibility  become  important.  Branched  compounds  usually  have 

much  lower  values  of  E,  than  nonbranched  compounds  of  the 

d  r 
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same  size. 

e.  Nature  of  the  Polymer.  From  the  hole  theory  of 

diffusion,  the  rate  of  diffusion  depends  on  (a)  the  number 

and  size  distribution  of  pre-existing  holes,  and  (b)  the  ease 

of  hole  formation.  It  is  physically  obvious  that  increased 

crystallinity  will  reduce  D.  Also,  if  cross  linking  increases 

the  diffusion  constant  decreases  and  the  activation  energy  of 

diffusion  will  increase.  It  is  also  relevant  to  mention  that 

bomb  of  these  processes  would  be  expected  to  decrease  the 

configurational  entropy  of  the  polymer  and  that  the 

correlated  decreases  in  D  would  be  qualitatively  consistent 

IS 

with  Barker’s  entropy  correlation  theory,  according  to  which 
ar.y  structural  change  which  reduces  the  configurational  entrap 
3„  of  a  complex  molecular  system  (the  polymer)  will  lead  to 

* 

an  approximately  equal  increase  in  the  activation  entropy  '3 
for  structurally  controlled  rate  processes.  The  addition  of 
plasticizer  to  a  polymer  decreases  the  cohesive  forces  between 
the  chains,  resulting  in  an  increase  in  segmental  mobility. 

It  is  clear  that  the  rate  of  diffusion  will  increase  when 


plasticizer  is  present. 


f.  Glass  Transition  Temperature.  Most  of  the  polymers 
show  a  sharp  change  of  Arrhenius  slope  at  their  glass  transi¬ 
tion  temperatures.  Above  th>.  gioou  t  oo 7  ore  tore  there  is  a 
greatly  increased  segmental  nobility  allowing  for  a  larger 

zone  of  activation,  leading  to  the  observed  increases  in  both 

19 

energy  and  entropy  of  activation.  Often,  below  the  glass 
temperature  diffusion  becomes  anomalous.  Sorption  curves  may 
shew  Non-Fickian  character.  D  may  show  time  and  thickness 
dependence . 

g.  Fickian  Characteristics.  The  mathematical  solution 
mentioned  in  section  I.3.2.C.  (Sorption  and  Desorption  Method) 
in  this  thesis  is  based  on  Fick's  law.  Sorption  curves  having 
characteristics  expected  from  the  specified  solutions  are 
customarily  called  Fickian  (or  normal)  type  sorption.  The 
following  is  a  summary  of  important  Fickian  sorption  features 
that  have  been  made  available  by  the  mathematical  studies  of 
Crank  and  associates. 

(a)  Both  absorption  and  desorption  curves  are  linear  in 
the  initial  stage.  For  absorption,  the  linear  region  extends 
to  over  60&  or  more  of  M  ,  where  (M  /  M  )  is  the  amount  of 

00  W  o 

vapor  absorbed  per  gram  of  dry  polymer  until  the  sorption 
equilibrium  is  reached.  For  D(C)  increasing  with  C  the 
absorption  curve  is  linear  almost  up  to  the  final  sorption 
■equilibrium. 


(b)  Above  the  linear  portions  both  absorption  and 
desorption  curves  are  concave  coward  the  abscissa  axis, 


irrespective  of  the  form  of  D(C). 

(c)  When  the  initial  concentration  and  the  final 
concentration  are  fixed,  a  series  of  absorption  curves 
for  films  of  different  thicknesses  are  superposable  to  a 
single  curve  if  each  curve  is  replotted  in  the  form  of  a  re- 
ducod  curve,  i.e.  (M  /MQ)  is  plotted  against  t2/b.  This  also 
applies  to  the  corresponding  series  of  desorption  curves. 

(d)  The  reduced  absorption  curve  so  obtained  always 
lies  above  the  corresponding  reduced  desorption  curve  if  D 
is  an  increasing  function  of  C  between  Ct  and  C^.  Both 
reduced  curves  coincide  over  the  entire  range  of  t  when  D  is 
constant  in  this  concentration  interval.  The  divergence  of 
the  two  curves  becomes  more  marked  as  D  increases  more 
strongly  with  C  in  the  concentration  range  considered. 

(e)  For  absorptions  from  a  fixed  to  different  C^'s, 

the  initial  slope  of  the  reduced  curve  becomes  larger  as  the 

concentration  increment  C_-C^  becomes  larger,  provided  that 

D  increases  monotonoca 1 ly  with  C  in  the  range  considered. 

This  also  applies  to  the  reduced  desorption  curves  which 

start  from  different  C-'s  to  a  fixed  C-. 

i  r 
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Chapter 


THEORIES 


In  this  chapter  a  number  of  theories,  which  have  been 
considered  in  the  analysis  of  the  experiments,  are  outlined 
Although  each  of  the  theories  discussed  provides  some 
useful  qualitative  insights,  a  self-consistent  combination 
of  the  theories  (mainly  the  DiBenedetto-Paul  theory  and  the 
dual  sorption  theory)  has  been  adopted  for  quantitative 
use  with  the  data.  It  also  should  be  mentioned  that  there 
is  considerable  overlap  of  the  concepts  in  the  various 
theoretical  models. 


Molecular  and  Statistical  Mechanical  Models- 
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1.  Molecular  Model.  DiBenedetto  and  Paul"  assumed  the 


sorbed  gas  in  an  ordered  region  of  a  polymer  is  contained 
within  a  cell  composed  of  four  parallel  segments  containing 
n  structural  units  per  segment.  The  activation  energy  for 
diffusion  is  related  to  the  variation  in  the  average  poten¬ 


tial  energy  of  a  molecular  interaction  which  accompanies  a 


volume  chancre  of  An(ir/4)r  ,  where  r  is  the  diameter  of 

'  g  g 


diffusant  and  \  is  defined  in  Fig.  2.  First  they  applied 
the  Leonard-Jones  6-12 (Mie)  potential  function 


=  c  [  ( o  /  o  ’ ) 


n  1  \  ^  ^  _ 


2  (  o* / o ’ ) 6 ] 


(20' 


to  this  system  using  a  model  in  which  a  single  center 


’  >0 


infracts  with  its  neighboring  chain  of  21  canters 
(arbitrary)  to  give  the  total  interaction  energy  <p  as 
below 

<b=£*  y  [  ( p  *  /  p !  )  1 2  -  2(p*/p:)S]  (30) 

c  i=-10  1 

*  ★ 

where  z  and  p  are  energy  and  distance  parameters 
respectively;  and  o'  is  the  distance  between  centers. 

Figure  (2)  gives 


2  2^2 
i>l  =  Pq  +  1  (a)  /  i 


(0  ,  ±1 , ±2 


) 


(31) 


where  ,\  is  again  the  chain  length  per  center  (twice  the  value 
of  A  in  ref.  20)  and  p  is  the  nearest  neighbor  distance. 
Substitution  into  equation  (30)  gives 


bc  =  £  LC  (p  /P^)^  -  2D(p  /Pq)  '3 


where 


C  = 


1  J-  ^  2  2-0 

1+2  )  [1  +  i  (A /o')  ] 

i=l  ° 


i=  10 


.  2 


~  i  r  2  ll 1  +  i  ( A/p ' ) “] 
i=l  u 


Now 


t  *  k  i  i  /  o  *  5/2 

?,  =  (  -  ■  /  A  )  [  0 . 7  7  ( v  /v)^'"  -  2.32  (v  /v)  j 


(32) 


(33 


’(.'')“  is  the  effective  volume  per  center,  which 


where  v 


came  from  the  assumption  that  each  center  is  constrained  to 

*  *  2 

move  in  plane  perpendicular  to  the  chain  axis  and  v  =  ,\p 

Diffusion  is  assumed  to  occur  when  an  extra  volume  of 
2 

An(7r/4)r  added  to  the  normal  effective  volume  nv. 

g 

2 

nv  =  Xn(i/4)r  +  nv.  (34) 

a  g 

The  activation  energy  is  given  by 


AEd  =  >h  (A$) 


with  A(p  equal  to  four  times  the  change  in  Eq.  33  correspond¬ 
ing  to  the  volume  change.  Equation  (33)  is  multiplied  by 
4  because  of  the  interaction  between  four  segments  of  this 
hypothetical  cell.  The  final  result  of  DiBenedetto  and  Paul 


★  ★ 

AEd  =  (4nNAc  IX)  •  p  { 0 . 77  [  (— )  1X/ -  -  t^-)  '  ] 


,v\s /: 


-  2.32[<— )3/z  -  (^)3/2]  } 
v  v 

a 


2.  Statistical  Mechanical  Model .  Pace  and  Datyner 


21,22 


used  the  same  assumption  as  DiBenedetto  and  Paul.  See 
Figure  (3).  An  average  chain  co-ordination  number  of  4  is 
assumed  rather  than  6  as  for  crystalline  closest  packing  of 
rods.  The  diffusant  molecules  may  move  either  (a)  by 


sliding  longitudinally  down  the  x-axis  parallel  to  the 


interstitial  "channels"  of  a  chain  bundle  or  (b)  by  "jumping" 
at.  right  angles  to  this  direction  when  adjacent  chains  are 
S'.: ; iiciently  separated.  One  may  readily  show  that  process 
(aj  is  generally  much  faster  than  process  {by  and  so  the 
penetrant  proceeds  until  a  barrier  {entanalment ,  crystallite, 
etc.)  is  reached  whose  penetration  requires  more  energy 
than  chain  separation.  Process  (b)  is  then  rats  determining 
and  gives  rise  to  the  observed  activation  energy  of  diffusion 
and  jump  frequency.  Process  (a)  establishes  the  root  mean 
square  displacement  (L)  between  jumps. 

Pace  and  Datyner  used  a  modified  form  of  the  equation 
from  DiBenedetto 

f(z)  =  (e*p*/A2) (0.77 (p*/z) 11  -  2.32(p*/z)5)  (37) 

This  equation  is  like  equation  (33)  except  for  an  extra 
division  by  A,  the  spacing  between  attractive  centers,  so 
that  f(z)  is  the  potential  energy  per  unit  length  of  a  chain. 

Neglecting  the  possibility  of  segmental  "flips"  and 
assuming  that  chain  bending  is  distributed  over  a  number  of 
backbone  bonds,  one  may  show  that  the  minimum  energy,  AE, 
necessary  to  produce  a  separation,  d,  between  adjacent  chains 
(neglecting  any  other  elastic  distortion  of  the  environment) 
is  given  to  good  approximation  by 

oo 

AE  =  min  j  (f(z)  +  ( 3/2 ) (d2z /dx" ) - - f ( o ) dx 


(33) 


subject  to 


■k 

z  =  p+  d  at  x  =  0 

z  =  o  at  x  =  ±<=°  (39) 

The  second  term  in  the  integrand  gives  the  intramolecular 
(chain-bending)  energy  and  x  measures  distance  along  the 
chains.  B  is  an  effective  chain-bending  module  and  is 
estimated  from  rhe  polymer  backbone  geometry  and  the  shape 
of  the  bend  rotation  potential  minima.  By  a  suitable 
linearizing  approximation,  one  may  solve  equation  (38)  to 
give 


E  = 


5. 23  (-£-)  1/4  (£-§-)  3/4  (0. 077  [  (£-)11(p 

u  ,  z  p 

A 

*  * 

*  n  1  0  n  c: 

-  P  ( — *  - )  1  -  0.58 [ {^— )  ( 0  -  4d ' ) 

p  +  d  p 

*,  o*  i  4,  t  3  /  4 

-  p  (— - )  i  ) 

o  +  d 


-  lOd ' ) 


(40) 


d'  =  d  +  p  -  p 


(41) 


Pace  and  Datyr.er  assume,  both  above  and  below  T  , 

g 

that  thermal  expansion  occurs  by  a  uniform  average  increase 

in  interchain  spacing.  Therefore,  neglecting  expansion 

* 

along  the  chain  one  has  (since  p  /p  C  1.0) 


i 


(d/dT) (p  /p) 


m  is  the  mass  per  L-J  center  (daitons) ,  length  units  ar 
nr.;,  and  energy  units  J/mol,  d  eeuais  d  (taken  from  gas 
viscosity  measurements) . 


For  polymer  chains,  such  as  those  indicated  in  Fig. 

*  * 

bulky  side  groups,  d  is  given  by  d  =  d  -  (p  -  Pc) • 


Fig.  4  Diagramatic  view  of  two  hypothetical  closely 

packed  helical  polymer  chains  that  possess  closely 

spaced,  bulky  side  groups.  The  solid  cylinders  define 

the  effective  outer  surface  of  the  chain  "cores " ;  the 

dashed  cylinders  are  defined  by  the  average  Leonard- 

Jones  parameter P*.  There  is  considerable  free  space 

’  22 

within  the  structure. 


29 

Equation  (40)  is  the  minimum  energy  to  separate  two  chains  by 
an  amount  d,  as  determined  by  their  stiffness  and  mutual 
attraction,  and  neglects  any  displacement  of  compression  of 
the  environment.  It  is  therefore  a  lower  bound  for  the  true 

activation  energy.  The  expression  contains  no  adjustable 

-2 

parameters  but  L  in  equation  (44)  is  net  predictable  within 
the  limits  of  the  present  theory. 

B.  Free  Volume  Theory 

2"* 

Cohen  and  Turnbull  “  treated  diffusion  as  translation 
of  a  molecule  across  the  void  within  its  cage.  In  other 
words,  diffusion  is  a  result  of  redistribution  of  the  free 
volume  within  the  polymer  and  to  some  extent  the  diffusar.f 
is  carried  along  as  the  cage  moves. 

The  free  volume  of  a  given  diffusant  molecule  is  de¬ 
fined  as  the  volume  within  its  cage  less  the  volume  of  the 
molecule.  Let  the  total  free  volume  Vc  be  divided  into  i 

i. 

small  regions  or  packets;  each  molecule  in  the  ith  region 
has  an  average  free  volume  v^.  Nn.  is  the  number  of  diffusan 
molecules  having  free  volume  in  the  ith  region. 

In  order  to  correct  for  overlap  of  free  volume,  a 
numerical  factor  m  is  introduced  to  give 

m  7  N  .  v  .  =  V  - 
.  ii  t 
l 


(46) 


The  value  of  m  should  lie  between  h  and  1. 


The  average 


free  volume  is  given  by  v^  =  V^/M  and  the  distribution 
function  P  ( v 1 )  is 

P(v')  =  (tn/vf}axp(-mv'/vf)  .  (4 

•k 

The  probability  of  finding  a  hole  of  volume  v  or  larger 
s 

CO 

? ( v  )  =  J  ? ( v ’ ) d v '  -  exp { -mv  /v ^ )  (4 

v  J" 

Cohen  and  Turnbull  assumed  a  molecule  moved  trans¬ 
lational  ly  across  the  void  within  its  cage  with  gas 
kinetic  velocity  u.  The  distribution  of  this  molecule  to 
the  diffusion  coefficient  is 

k 

D  ( v )  =  ga(v'Ju  for  v'  >  v  (4 


where  a(v')  is  roughly  the  diameter  of  the  cage  and  g  is 


is 


D  =  j  D ( v ' ) P ( v ' ) dv '  =  ga  uexp(-mv  /vc) 


(5 


★ 

where  a  equals  approximately  the  molecular  diameter. 

It  is  chon  further  assumed  that  the  free  volume  vr 
be  expressed  bv  an  aooroximate  exoression 


3 


where  Tq  is  the  temperature  at  which  the  free  volume 

disappears  a  and  v  are  the  mean  values  of  the  coefficient 

m 

of  thermal  expansion  and  the  molecular  volume  over  the 
temperature  range  of  interest,  respectively,  substitution 
of  equation  (51)  into  equation  (50)  gives 


*  *  _ 

D  =  ga  uexp(-mv  / av^ ( T  -  TQ) ) 


C.  Earner's  Zone  Theorv 


Barrel 


based  his  theorv  on  the  idea  that  a  diffusina 


molecule  moves  to  successive  positions  as  segments  of  the 
surrounding  polymer  chains  fluctuate  in  an  activated  zone 
within  which  the  activation  energy  is  distributed  through 
many  degrees  of  freedom.  The  expression  for  the  diffusion 


coefficient  is 


f  'max 


D  =  ( v  /  2 )  ,\  '  l  P  [(E/RT)-  (1/  (f  ’-1)  !  }j  exo(-E/RT$  (5 

f=l  r 

where  v  is  the  thermal  vibration  frequency  of  the  penetrant 
molecule.  A'  is  the  jump  distance,  f'  is  the  number  of 
degrees  of  freedom  involved  in  a  particular  movement,  f'max 
is  the  value  of  f'  for  which  the  expression  in  the  square 
bracket  has  a  maximum,  E  is  the  total  energy  of  1  mole  of 
the  activated  zones  under  consideration,  and  Pf  is  the 
probability  that  the  £'  degrees  of  freedom  will  co-operate 
in  a  diffusion  step.  The  complete  expression  in  the  square 


-  ;•  /■,  v>  /  w v  -v/i  * •. i 


•v  *>, 


brackets  is  the  probability  that  the  energy  E  is  distributed 
over  £'  degrees  of  freedom. 


D.  Transition  State  Theory  (Absolute  Race  Theory) 

25 

Eyring  proposed  a  model  that  the  molecule  vibrating 
within  an  energy  well  had  a  tendency  occasionally  to  jump 
over  an  energy  barrier  to  the  new  energy  well,  causing 
diffusion  to  occur. 

The  expression  for  D  is  given  by 


D  = 


gA'2UT/h)e4S‘/V5H*/RT 


(:J) 


J 

* 


} 

■i 


where  AS*  and  AH*  are  activation  entropy  and  activation 
enthalpy  respectively.  A'  is  jump  distance,  k  the  Boltzmann 
constant,  h  Planck's  constant,  and  a  represents  a  geometri¬ 
cal  parameter  which  reflects  the  dimensionality  of  the 
diffusion  process.  For  isotropic  diffusion  in  three 


dimensions,  g  =  1/6. 

The  rate  of  jumping  (kT/h)e  e  is  given 

by  the  product  of  the  attempt  frequence  kT/h  and  the 

probability  per  attempt  of  hopping  over  the  potential 
-  A  G*/RT 

Darner,  e 

In  more  recent  variants  of  the  rate  theory,  the  jumping 
probability  is  enhanced  by  quantum  tunneling  through  the 


barriers . 


' '■  .••V-  t'  .  'lSV»\v .-•VyV. *.'-V  V  .'  A  ->-■  y- 


Dual  Sorption  Theory 


25 


Wolf  R.  Vieth  and  Mary  A.  Amini  formulated  the  "dual 
sorption  theory"  to  explain  certain  classes  of  negative  and 
positive  deviations  from  Henry's  law  which  are  frequently 
observed  in  sorption  plots  of  penetrants  in  polymers.  The 
dual  modes  are  either  that  the  penetrant  molecule  is 
normally  dissolved  and  free  to  diffuse  or  that  it  is 
immobilized,  as  in  a  sink  or  well.  It  is  the  second  prcces 
which  gives  rise  to  deviations  from  normal  behavior. 


1.  Negative  Deviation  From  Henry's  Law.  The  negative 
deviations  can  be  associated  with  site  binding  and  low 
penetrant  levels,  where  some  of  the  penetrant  molecules  are 
immobilized  on  sites  or  in  sinks  in  the  matrix.  This  would 
be  characteristic  of  a  glassy  polymer. 

The  over  all  concentration  of  sorbed  vapor  can  be 
written  as  the  sum  of  two  terms: 


C  =  Cd  +  Sh  (55) 

where  represents  ordinary  dissolution  and  S  represents 
sorption  in  the  microcavities  or  holes  (S^) ,  Vieth  and 
Amini  explained  the  sorption  behavior  by  using  a  modal 
based  on  Henry's  law  and  the  Langmuir  isotherm, 

C^j  =  K^p  (Henry's  Law)  (56) 


(Langmuir  Isotherm) 


(57) 


S.  =  S.'bp/  (1  +  bo) 
n  n 

where  s,'  =  the  hole  saturation  constant 
h 

b  =  the  hole  affinity  constant 
At  lew  pressures,  where  bp  <<  1,  the  sorption  isotherm 
tends  to  linearity 


c  =  (Kd  +  biV)p  (58) 

In  the  limit  of  high  pressures,  where  bp  >>  1,  the  solu¬ 
bility  in  the  holes  reaches  the  saturation  limit,  . 
Therefore,  the  isotherm  should  once  again  become  linear 


C 


K.p  +  S' 
dr  n 


(59) 


Thus,  a  plot  of  C  versus  p  will  have  two  linear  regions 
with  a  connecting  nonlinear  region.  Vieth  and  Am ini 
applied  these  relations  to  diffusion.  To  analyze  the 
diffusion  process  for  dual  mode  sorption,  a  more  complex 
model  is  necessary.  A  differential  mass  balance  is  per¬ 
formed  for  penetrant  diffusing  across  a  fiat  film.  One 
term  (left  term  in  equation  60)  allows  fox-  the  rate  of 
accumulation  by  simple  dissolution;  the  other  (right  term  in 
equation  SO)  for  the  fact  that  the  film  is  also  entrapping 
an  amount  of  the  penetrant  and  holding  it.  The  result  is  a 
modification  of  Fick's  law. 


C-V-v  , 


(60) 


D-- 


3x ' 


at  h 


(SK  +  Ch 


For  sorption  of  a  gas  in  a  glassy  polymer  the 


■ ir 


urn  partial  pressure  of  the  gas  at  any  position  ir. 


tiie  film  is  the  same  for  both  types  of  sorption.  By 
equating  the  pressure  from  the  two  formulations,  a  relation 
ship  is  found  between  the  concentrations  of  the  bound  and 
the  free  species. 


Sh  -  [<S£b/Kd)Cd)/ll  ♦  !b/KdICdi 


(61) 


Substituting  this  expression  for  Sj,  into  equation  (60) 
yields  the  diffusion  equation  for  dual  mode  sorption. 


o 

3  C,  3C  . 
d— 4  =  [i 


shlb/!:d> 


dX 


3r 


(1  -  tfc/K.) 


(62) 


d'^d' 


At  low  pressures  equation  (61)  reduces  to  a  linear 


relationshio  between  S,  and  C-, 


u 


Sh  =  RCd 


(63) 


wnere 


R  =  bS.1  /K  - 


(6-4) 


of  Pick' 


lav;,  wn«. -re 
•> 


iusivity  is  defined. 


err. 


.  t 


( 6  -  ) 


where  D  ^ c  -  the  stfective  diffusivity 
erf  1 

=  D/(l  +  R)  (66) 

The  measured  diffusion  coefficient  is  less  chan  the 
true  diffusion  coefficient;  that  is,  the  attainment  of 
sorption  equilibrium  is  slowed  down  by  the  immobilization 
process . 

At  higher  pressure,  the  partial  differential  equation 
for  diffusion  is  nonlinear.  Numerical  solutions  can  be 
found  for  the  equation  by  using  finite  difference  techniques 

For  normal  sorption  processes  occurring  in  finite 
baths  of  penetrant,  the  dimensionless  pressure  decay,  is 
a  simple  function  of  the  square  root  of  dimensionless  time  6 

8  =  Dt/L2  (67) 

where  2L  =  the  thickness  of  membrane 


4>  =  <?j_  ~  P)/(?j_  “  Pf)  (68) 

v/here  =  the  initial  pressure 
P|r  =  the  final  pressure 

The  pressure  decay  is  plotted  versus  the  square  root  of 

'•  =  Ml  +  (S^b/Kj  /  (1  +  bp)2)'1.  (69) 


Vieth  and  Anum  round  three  numerical  solutions  by  varying 
the  sorption  parameters,  b,  5,'  ,  and  7.^  over  a  wide  range. 
Expressed  graphically  all  of  the  data  were  correlated  with i 
a  narrow  band.  Therefore,  the  mean  of  the  band  was  taken 
as  the  correlation. 

2.  Positive  Deviation  From  Henry's  Law.  In  this  case,  the 
polymer  network,  swells  to  expose  more  sites,  increasing  the 
sorption  level  synergistical Iv .  Now,  the  site  density  is 
not  constant  when  swelling  occurs;  it  changes  with  the 
extent  of  sorption.  A  clustering  phenomenon  displays  a 
significant  deviation  from  Henry's  law,  which  predicts  no 
cluster  formulation  at  ail.  A  more  detailed  discussion  of 
this  point  is  not  appropriate  in  this  thesis. 

The  result  is 


Deff  =  i^/U  +  Ik)]  ( 1  +  ( 3lny/31nC)  ) 


(70) 


In  this  case,  the  true  diffusion  coefficient  is 
modified  by  two  factors;  one  for  the  second  mode  of  sorpt^o 
which  reduces  it,  and  another  for  network  swelling,  which 
increases  it. 


Barker's  Entroov  Correlation  Thoorv 


18 


The  central  point  ot  the  thoorv  is  that  in  a  struc¬ 
turally  complex  system  undergoing  a  thermally  activated 


v 


-r.  w.  • 


v£" 


5 

rate  process,  any  organizing  influence  which  produces  a 

reduction  -AS  in  the  configurational  entroov  will  be 
c  ^ 

★ 

accompanied  by  a  correlated  increase  A (AS  )  in  the  activati 
encropv  associated  with  the  rate  process.  Thus, 


A  (AS  )  =  -AS 

C 

*  * 

AS  - .  .=  AS AS  . 

rinal  initial  c 


(71) 


(72! 


The  experimental  support  for  this  theory  was  obtained 
by  studying  permeation  of  various  gases  in  a  series  of 
oriented  polyalkyl  methacrylates.  The  conclusions  drawn 
were  that  in  a  ncr.crystalline  polymer,  when  the  configura¬ 
tional  entropy  Sc  is  reduced  by  an  ordering  process  such  as 
stretching,  an  additional  amount  of  local  disorder  must 
occur  in  the  neighborhood  of  a  diffusing  molecule  before  a 
unit  diffusionai  step  can  occur,  so  that  the  decrease  in 
the  configurational  entropy  is  equal  to  the  increase  in 
the  activation  entropy  for  diffusion.  It  may  be  seen  that 
this  theory  is  in  at  least  qualitative  agreement  with  many 
ether  facts.  For  example,  crystallization  reduces  the 
configurational  entropy  so  that  one  would  expect  a  large 
increase  in  the  activation  entropy  and  also  a  large  increas- 
in  the  activation  energy.  The  competition  between  those 


two  effects  appears  to  be  such  that  without  the  entropy 


39 

correlation,  the  value  of  D  would  be  much  more  strongly 
reduced  than  it  actually  is. 

"k 

C.  Ionic  Conduct ion  in  Polymers 

In  a  given  polymer  system,  it  should  be  considered 
that  electrical  conduction  may  be  due  to  the  movement  of 
any  of  the  following  types  of  charge  carriers:  electrons, 
holes,  protons,  bare  elemental  ions,  molecular  ions,  and 
solvated  or  hydrated  ions.  From  the  viewpoint  of  diffusion 
experiments,  ionic  conduction  should  play  the  most  important 
role  in  this  sort  of  study.  Barker's  "Local  Structure 
hypothesis"  provides  a  useful  concept  to  approach  the 
problem.  The  idea  is  that  the  actual  concentration  of  ions 
needed  to  provide  a  typical  level  of  conductivity  (a)  is 
suite  small  compared  to  the  number  of  mers  present  in  the 
same  volume  and  furthermore,  that  she  electric  field  near  a 
given  ion  is  very  large,  so  that  a  measurement  of  a  proves 
a  small  volume  fraction  of  the  sample  in  regions  that  may 
be  atypical  due  to  one  presence  of  the  ions. 

The  ionic  conductivity-  of  polymers  can  originate  either 
from  ionically  dissociable  groups  in  the  pure  polymer  or 
from  the  dissociation  of  an  ionic  impurity  in  the  polymer. 
The  latter  case  is  much  more  common.  Water  plays  an 

★ 

The  close  collaboration  of  D.  Y.  Chen  in  this  section 
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important  role  in  conduction  in  polymers  according  to 

2  8 

Barker  and  Sharbaugh's  "Weak  Electrolyte  Model."  The 
enlacing  ionic  impurity  will  be  only  partially  dissociated 
as  a  weak  electrolyte  into  aniens  and  cations  due  to  the 
presence  of  water.  The  pure  water  itself  will  also  provide 
H+  and  OH  (or  H^O’)  ions.  These  dissociated  ionic  species 
will  move  toward  the  electrodes  when  an  external  -lectric 
field  is  applied.  The  process  gives  an  electrical  con¬ 
duction  current. 

One  approach  used  to  deal  with  the  problem  of  electri¬ 
cal  current  due  to  ionic  diffusion  is  the  concept  of  an 
activation  barrier. “  For  a  simple  ionic  solid  composed 
of  atomic  ions,  the  diffusing  ion  would  move  either  by 
exchange  with  adjacent  vacancy  or  by  being  excited  to  an 
interstitial  position  from  which  activated  diffusion  become 
easier.  For  a  general  diffusion  mechanism,  the  number  per 

unit  volume  of  ions  (n)  in  sites  from  whicn  motion  is 

-W/kT 

possible  is  n  =  Me  ,  where  N  is  the  total  number  per 

unit  volume  of  ions,  W  is  the  formation  energy  of  sites.  I 

the  activation  carriers  for  jump  itself  is  U  and  the 

frequency  of  vibration  is  v,  the  number  of  jumps  per  second 

for  each  of  n  potential  diffusive  centers  in  the  absence  of 

“  (J  /  k  T 

an  applied  field  is  ve  ^  .  Thus,  the  total  jumps  per 

— ( w+U ) /kT 

second  is  JQ  =  Nve  '  A 1 g .  The  effect  of  an  applied 


field  E  is  to  raise  the  activation  barrier  in  the  direction 
of  the  field  by  eEA'/l,  where  X'  is  the  average  diffusion 
jump  distance,  and  to  lower  the  barrier  in  the  opposite 
direction  by  the  same  amount  (Fig.  5) .  The  factor  g  is  a 
term  involving  geometrical  and  entropic  features  of  the 
system.  The  net  number  of  jumps  per  second  in  the  directior 
opposed  the  field  is 

J  =  J_  -  J+ 

=  \'gNve_W/kT(exp(- (U-eEX'/2) /kT)-exp(- (u+eEX '12) /kT) 

=  2a 'gNve- (W+U> /kTsinh (eEX’/2kT) 

=  2 JQsinh (eEA ' / 2kT)  (73) 

The  current  density  (j)  due  to  this  diffusion  process 
is  j  =  Jq  =  2NvqgX ’e"  (W+U) /xTsinh  (eE.\ '/2kT)  (74) 

where  q  is  the  charge  per  ion. 

The  implications  from  the  above  equation  are: 

(i)  the  conduction  is  non-ohmic  at  high  values  of  the 
applied  field 

(ii)  j  increases  with  temperature  T. 


CaaDter  III  INSTRUMENTATION  AND  EXPERIMENTAL 


PROCEDURE 


A.  X-Ray  Energy  Dispersive  Analysis 


1 .  Instrumentation .  X-ray  analysis  has  become  widely 
accepted  as  a  powerful  tool  for  qualitative  and  quantita¬ 
tive  analysis  of  the  distribution  of  chemical  elements 
because  it  is  nondestructive,  is  direct  in  sample  preparation 
techniques,  can  analyze  nearly  all  of  the  elements  of  the 
periodic  table,  and  is  very  accurate.  It  is  unique  among 
analytical  instrumental  methods  in  that  it  provides  a  con¬ 
siderable  amount  of  information  quickly.  Within  minutes, 
useful  data  concerning  the  atomic  composition  of  the  sample 
for  the  elements  with  atomic  numbers  from  Na  up  can  be 
generated . 

a.  General  Principles.  When  an  X-ray  photon  or 
electron  beam  from  a  source  interacts  with  the  inner  shells 
of  an  atom,  an  electron  moves  from  one  shell  to  a  higher, 
but  unstable  level.  When  this  atom  returns  to  its 
stable  position,  an  X-ray  is  released  that  is  characteristic 
of  the  particular  element.  These  X-rays  are  sensed  by  a 
Si  (Li)  detector,  amplified,  digitized,  and  stored  in  memory 
f-'.r  use  by  a  microprocessor.  Such  data  are  then  reduced  to 
answers  that  are  meaninaful  to  the  occrator. 
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b.  Spectrometer .  There  are  two  types  of  X-ray 

spectrometers.  One  is  wavelength-dispersive,  and  the  other 

one  is  energy-dispersive.  The  energy  dispersive  X-ray  (EDX) 

systems  are  lower  in  cost,  but  may  not  provide  sufficient 

resolution  or  may  not  have  the  ability  to  handle  high  con- 
30  31 

centrations.  ’  ^  The  wavelength-dispersive  X-ray  (WDX) 
spectrometer  has  better  resolution,  but  it  needs  time  to 
scan  the  entire  spectrum  of  interest  with  the  diffraction 
crystal.  The  systems  analysis  time- typically  between  50  s 
and  3  min  for  a  single  sample  with  as  many  as  eight  measured 
elements  is  somewhat  longer  than  more  complicated  instru¬ 
ments.  Also  this  kind  of  spectrometer  will  be  affected  by 
the  geometry  of  the  sample.  So,  the  proper  alignment  be¬ 
tween  the  collimating  slits,  the  diffraction  crystal,  and 
the  direction  should  be  established.  The  spectrometer  used 
in  this  experiment  is  of  the  EDX  type.  Figure  (6)  is  a 

sc.nematic  diagram  of  the  energy  dispersive  X-ray 
32 

Spectrometer .  The  detector  is  a  semiconductor  whre.:  can 

distinguish  one  chemical  from  another  on  an  energy  selection 

basis.  The  typical  Si (Li)  spectrometer  is  diagramed  in 
3  3  3  4 

Figure  (7).  '  The  Si  (Li)  detector  is  essentially  a 

reverse-bias  diode  which  acts  as  a  transducer;  converting 
the  energy  of  the  X-ray  photons  to  electron-hole  pairs. 

The  electron-hole  pairs  are  swept  out  of  the  detector  by  the 
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applied  bias  voltage  and  the  total  charge  produced  in  the 
current  pulse  is  directly  proportional  to  the  energy  of  th< 
X-ray  that  was  absorbed  in  the  detecror.  The  Si (Li)  diode 
and  the  field  effect  transistor  (FET) ,  which  is  the  first 
amplifier  stage,  are  cryogenically  mounted  to  red1  '"•a  the 
electronic  noise  of  the  devices.  The  detector  has  good 
operational  performance  at  the  boiling  point  of  nitrogen, 
-L30°C,  but  the  FET  must  be  warmed  70°C  above  this  for 
optimal  performance,  necessitating  the  heater  resistor. 

The  amplifier  increases  the  signal  level  from  the  FET 
and  also  reduces  noise.  Before  the  signal  goes  to  the 
multichannel  analyzer,  it  enters  the  base  line  restorer 
which  maintains  a  base  line  between  pulses.  The  base  line 
restorer  output  is  passed  to  an  analog-to-digital  converte 
(ADC) .  Each  time  an  X-ray  pulse  is  processed  by  the  ADC  a 
number  is  generated  that  represents  the  channel  to  be  in¬ 
cremented  by  one,  thereby  resulting  in  an  intensity  versus 
energy  distribution.  A  magnetic  core  memory  accumulates 
the  number  of  times  each  energy  interval  is  detected  and 
the  information  is  displayed  on  a  cathode  ray  tube  (CRT) . 

c.  Applications.  The  applications  include  scanning 
electron  microscopes,  electron  microprobes,  transmission 
electron  microscopes  and  X-ray  fluorescence  analysis.  The 
one  used  in  this  experiment  is  the  scanning  electron 


microscope-  Figure  (8)  shows  the  interface  between  the 
scanning  electron  microscope  and  energy  dispersive 
spectrometers.  The  SEM  has  a  very  small  electron  beam. 

For  this  reason  the  regions  of  the  sample  examined  can  be  as 
small  as  1-2  pm  in  diameter. 

d.  System  Performance.  Because  the  width  of  the  band 
pass  filter  and  the  fact  that  only  one  X-ray  can  be 
processed  at  a  time,  the  output  count  rate  is  less  than  the 
input  count  rate.  To  correct  for  this  difference,  the 
system  uses  a  special  livetime  clock  that  runs  only  during 
the  time  the  amplifier  is  able  to  accept  another  pulse  for 
processing.  When  an  X-ray  is  detected,  the  amplifier  be¬ 
comes  busy,  and  the  livetime  clock  is  turned  off.  The 
system  is  dead  during  this  period,  and  so  this  processing 
tame  is  called  deadtime.  Realtime,  the  normal  time  measure¬ 
ment,  is  the  sum  of  the  live time  plus  the  deadtime,  if  the 
deadtime  loss  is  over  20%,  the  livetime  correction  is  no 
longer  able  to  work  accurately. 

The  specified  resolution  for  the  system  was  153  eV  full 
width  at  half  maximum  (FWIIM)  at  5.89  eV.  The  reason  for 
using  this  is  that  the  theoretical  shape  of  the  curve  can 
be  predicted.  Usually  the  faster  the  amplifier  can  process 
counts,  the  poorer  its  resolution. 


Schematic  Showing  the  Interface  Between  an  Energy 
Dispersive  Spectrometer  and  a  Scanning  Electron  Micro 


2.  Experimental  Procedure .  Many  methods  ware  considered  as 
a  possible  means  of  introducing  diffusant  into  the  sample. 
Most  of  these  failed  because  of  the  severe  difficulty  in 
dealing  with  such  small  thin  films  (10  am- 4 0  qn  in  thickness? 
and  fine  fibers  (10  pm-40  urn  in  diameter) .  One  of  the 
methods  triad  is  shown  (Fig.  9. a.). 

Saits  such  as  NaCl,  Fe(NO^)^,.  .  .etc.  were  dissolved 
in  water  to  act  as  diffusants.  The  advantage  of  this 
arrangement  is  that  the  surface  tension  of  the  solution  and 
the  balance  pressures  (inside  tube  and  outside  tube)  create 
a  perpendicular  surface  to  introduce  diffusant  into  the 
fiber.  Furthermore,  the  evaporation  of  water  is  insig¬ 
nificant.  The  disadvantage  is  that  the  experiment  is  still 
complicated  by  surface  diffusion.  After  a  certain  time  of 
diffusing,  the  sample  was  glued  on  the  carbon  mount  using 
carbon  paste  (colloidal  graphite) .  The  reason  for  using 
carbon  was  to  avoid  peak  overlaps  between  the  specimen  and 
specimen  holder.  The  electron  microscope  only  allows 
conductive  specimens  to  be  detected.  Therefore,  the  specimen 

O 

was  coated  with  carbon  about  100  A  thick.  In  this  procedure, 
the  sample  was  exposed  in  a  vacuum.  We  assume  thac  the 
water  vaporized  in  this  step  and  that  diffusion  was  quenched. 
The  first  essential  problem  here  is  that  during  the  coating 
procedure  the  sample  was  heated  up  to  about  100°C  for  10 


Fig.  (9. a)  Method  of  introducing  diffusant 
into  PPBT  fiber. 
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Fig.  ( 9 . b )  Computer  output  of  qualitative 
analysis  method. 


minutes  for  the  process  to  finish.  Another  problem  is  that 
the  specimen  is  exposed  to  a  large  electric  field.  Con¬ 
sequently,  this  step  will  cause  a  severe  error  in  determining 
the  diffusion  rate  of  mobile  diffusant. 

The  coated  fiber  was  scanned  by  the  electron  beam  from 
an  edge  along  the  axis  at  each  distance.  If  diffusion  did 
occur,  the  excited  X-ray  could  be  analyzed  to  determine  the 
diffusion  rate. 

PP3T  fiber  25022-42-1  was  examined  by  EDAX  using  the 
qualitative  identification  program.  The  computer  output  is 
shown  (Fig.  9.b. ) . 

This  fiber  was  made  from  polyphcsphericacid  (P? A) 
solution.  The  residual  acid  corresponds  to  the  peak  of  the 
element  P.  The  peaks  of  Cu  and  Zr.  are  due  to  the  chamber 
materials . 


One  ?e(N0^)7  doped  PPET  F 
place  six  times.  me  results 
(10)  gives  the  computer  output 
analysis  method. 

The  proposed  explanations 
cause  errors  in  the  take  off  a 
inhomogeneous  structure  inside 


bar  was  examined  at  the  same 
hew  a  250i  deviation.  Figure 
of  the  semiquantitative 

are:  1.  the  rough  surfaces 

glo  (35°  here).  2.  The 
PPET.  3.  The  sensitivity  of 


EDAX  is  not  high  enough  tor  the  available  samples. 


to  the  change  in  sample  weight.  The  beam  is  always  in 
dynamic  equilibrium,  with  the  sum  of  the  moments  on  it  equal 
to  zero.  The  restoring  force  is  powerful  and  fast  so  than 
the  beam  appears  visually  to  be  locked  in  place.  3y  means 
of  accurate  potentiometers,  an  accurately  calibrated  voltage 
is  subtracted  from  the  voltage  across  the  coil.  A  dial  on 
the  potentiometer  reads  directly  in  milligrams,  correspond¬ 
ing  to  the  amount  of  voltage  being  subtracted  in  the  circuit 
The  access  of  coil  voltage  over  reference  voltage  is  then 
available  for  the  recorder. 

In  order  to  return  the  beam  to  exactly  the-  original 
position  after  a  weight  change,  a  slight  offset  of  the  beam 
is  required  to  cause  the  added  current  in  the  coil.  The 
racio  of  the  offset  which  would  occur  without  feedback  to 
that  which  occurs  with  the  servo  operating  is  called  the 
servo  loop  gain,  G.  It  is  a  dimensionless  quantity  in  exces 
of  1000.  This  factor  reduces  the  power  supply  changes  and 
speed  up  response  time . 

Figure  (14)  and  (15)  show  the  circuits  of  the  servo 
anoiifier  and  control  units  respectively.  The  force  induced 
by  the  phototube  current  is  limited  by  the  slit  opening. 
Optimum  adjustment  of  the  SLIT  potentiometers  L7  occurs 
when  the  maximum  up  and  down  forces  are  equal.  Due  to  this 
feedback  system,  drastic  changes  in  tube  characteristics 
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with  time  will  have  no  effect  on  balance  performance  within 
its  range.  When  the  beam  is  in  balance,  the  current  is  zero. 
Weight  added  to  the  sample  pan  will  cause  a  positive  current. 
On  the  other  hand,  removing  weight  from  the  sample  pan  will 
cause  a  negative  current.  The  maximum  current  is  ±3mA, 
maximum  voltage  about  ±9  volts,  correspondingly,  the  refer¬ 
ence  voltage  supplies  across  the  mass  potentiometers  P.29  is 
made  +  and  -  by  means  of  resistors  R30  and  R31,  so  that  the 
mass-reference-voltage  can  equal  coil-voltage  over  the 
useful  range.  Zero  net  weight  corresponds  to  0.5  on  the 
dial,  -  maximum  output  corresponds  to  0,  and  +  maximum  to 
1.0.  The  voltage  output  at  a  dial  seating  of  0.5  is  changed 
by  means  of  R32  to  compensate  for  variations  in  pan  and 
stirrup  weights,  etc;  it  is  basically  a  zero  control,  and  is 
called  "set  5,"  because  that  is  what  you  do  with  it.  The 
voltage  at  a  given  weight  is  adjusted  by  means  of  R22,  in 
series  with  the  torque  motor  coil.  The  more  resistance, 
for  the  same  current,  the  more  voltage.  This  is  similar  to 
the  Calibrate  controls  on  other  Electrobalances,  but  is 
called  SET  0/10  here. 

The  excess  of  coil  voltage  over  reference  voltage  is 
applied  through  R33  to  a  step  voltage  divider,  used  to 
establish  the  recorded  range.  The  actual  voltage  per 
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milligram  is  determined  by  adjusting  P.2  3  (Calibrate  Re¬ 
corder)  on  any  recorder  range. 

To  reduce  the  stray  forces,  from  vibrations,  esc.,  the 
output  voltage  can  be  applied  to  a  low-pass  filter,  depend¬ 
ing  on  recorder  characteristics.  If  more  filtering  is 
necessary,  add  more  capacitance  in  parallel  with  the  one  in 
the  balance.  Additional  capacitors  should  be  of  Mylar,  and 
should  be  added  in  increments  cf  5  or  10  pF  to  be  signifi¬ 
cant. 

Some  other  signal  can  be  introduced  to  the  recorder, 
by  connecting  its  +  lead  to  x  positive  of  the  recorder 
range  switch,  and  its  -  lead  to  2,  and  turning  switch  to  :< . 

b.  System  Performance.  There  is  a  residual  mechanical 
zero  shift  as  a  function  of  temperature ,  even  at  0.5000.  At 
no  load  it  is  about  -0.5  pg,'°C;  that  is,  the  balance  reading 
will  go  down  that  much  as  the  temperature  increases.  The 
effect  of  load  is  in  the  opposite  direction,  so  that  some¬ 
where  between  the  half-and  full-rated-load  the  effective 
temperature  coefficient  will  be  sere.  The  electromagnetic 
system  is  temperature  compensated,  so  that  the  MASS  dial  and 
recorder  calibrations  are  independent  of  temperature.  This 
compensation  is  full-,  effective  for  a  range  of  about  10°C 


around  ambient 
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Minimum  recommended  recorder  ranges  are  indicated  by 
the  small  pointers  on  the  back  of  the  mass  dial  range  knob. 
They  are  chosen  so  that  the  minimum  increase  of  mass  dial 
seating  will  be  11  on  the  finest  recorder  range.  If  the 
mass  dial  is  not  moved,  the  stability  of  the  reference  power 
supply  for  the  mass  dial  will  p'ermit  readings  to  1%  on  the 
recorder  for  2  ranges  below  the  minimum  recommended ,  at  full 
scale  on  the  mass  dial  (0  or  10).  If  the  mass  dial  happens 
to  be  close  to  5,  even  finer  ranges  can  be  employed. 

The  instruction  manual  #3124-1  for  #200-1  Cahn  electro- 
balance  model  RG  has  the  details  of  fine  points- of  operation, 
general  applications,  and  operating  procedures  etc. 

The  information,  such  as  opera ticn  in  a  vacuum, 
thermomolecular  flew  and  aerodynamic  effects,  is  very  useful 
for  the  person  working  in  this  kind  cf  area. 

In  the  experiments  of  absorption  carried  out  here, 
high  precision  was  needed  so  that  the  long-term  zero  draft 
and  cumulative  zero  drift  became  serious  problems.  The 
balance  needs  18  hours  to  warm  up  and  longer  may  even  be 
natter.  In  many  installations  the  blank  will  show  a  daily 
24-hour  cycle,  superimposed  on  a  steady  drift,  even  after 
all  temperature  effects  are  eliminated.  The  daily  cycle 
may  be  as  large  as  10  micrograms  peak-to-peak  at  no  load, 
and  15  micro-grams  at  full  lead.  The  steady  drift  rate  may 
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be  as  large  as  -s-0 . 1  microgram/hour .  The  explanation  for 
this  is  that  it  may  be  due  to  level  shifts  in  the  building. 
Other  factors  will  cause  drift,  too.  Periodic  blasts  of 
cold  air  from  an  air  conditioner  may  cause  periodic  errors. 
Turning  room  lights  on  and  off  may  produce  a  shift  of  as 
much  as  1-2  micrograms  in  balance  indication,  and  people 
walking  in  front  of  it  may  shift  'it  one  or  two  tenths  of  a 
microgram.  The  supply  voltage  change  of  ±  15%  will  cause 
a  few  tenths  of  a  micrcgram  shift. 

2.  Experimental  Procedure.  PPBT  films  23555-25-6  were 
cut  to  3.5  cm  long.  The  weight  was  ~7  mg.  Each  sample  was 
put  in  one  small  specimen  container.  Then,  they  were  trans¬ 
ferred  to  the  desiccator  for  48  hours  drying.  The 
desiccator  was  set  up  as  in  Figure  (16). 

After  the  samples  had  been  dried,  the  vacuum  was  re¬ 
leased.  During  this  releasing  process,  the  air  passed 
through  a  long  column  containing  desiccant  before  entering 
the  desiccator.  So,  the  moisture  of  the  air  was  absorbed. 
Then,  those  dried  samples  were  placed  in  the  electrobalance 
for  measurement. 

During  measurements,  the  mass  dial  was  adjusted  to  the 
value  where  the  recorder  went  to  the  zero  position,  and  it 
was  locked  at  that  moment.  The  recorder  was  in  the  highest 
sensitivity  range.  (2  ranges  below  the  minimum  recommended.) 


The  full  scale  (100%)  position  or  recorder  was  aown 
1  mg  at  the  5  mV-sefc.  Drifting  and  noise  became  a  bi 
b 1 am  but  we  had  no  choice.  The  raisons  for  dr if tin. 
have  been  discussed  in  the  section  on  System  Perrormance . 
Noise  was  caused  from  the  electrical  circuits.  Obviously 
we  were  dealing  with  the  instrument  near  its  margin  of 
capability.  There  was  one  other  important  factor  which  he 
to  be  considered:  the  recorder  did  not  record  the  mass 
change  from  instant  that  the  diffusion  began.  Figure  (17) 
is  based  on  a  curve  traced  from  the  recorder.  Points  were 

chosen  on  the  curve  and  plotted  as  relative  mass  increase 

x. 

vs  t2  in  Figure  (18).  The  curve  is  clearly  non-iinear  m 
the  initial  region.  If  we  assume  the  recorder  has  a  time 
lac  t_,  then  the  curve  is  naturally  based  on  (M  -  M,  )/ 

°  j.  t  t0 

(M  -  M.  )  vs  (t  -  tn)^.  Hence  consider 

'  »  r  0 

U 


(M,  -  M  ) / (M  -  ) 

u  fc0  ^  "0 


r  '  M  f  V  \  — 

-  N  1  *  4.  /  -  ‘  / 


(Mfc  /M  J  ]  /  [1  -  (Mt^ 


M./M  =  1  -  (3/rr  )  )  (2m  t  1)  exp[-D(2m+l)  r  t/o  ] 

t  --O  u ,, 

m=0 

(See  equation  (14)) 


According  to  above  equations,  we  can  fit  the  experimental 


curve  bv  chirmincr  t,  and  D.  The  best  fit  is  shown  on 


1.23  s  10’9cm2/s 


i 
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Fa cure  (19),  where  =  0.1  man,  D  =  1.05  x  10  cm/sand 
M.  /M  =  0.08954395433. 

Curves  such  as  the  one  in  Figure  (19)  were  drawn  using 
cue  21  step  curve-drawing  program  of  a  Hewlett-Packard  Desk 
Top  Microcomputer  with  a  7225A  plotter. 

C.  Sorption  Method  with  n  Sample  Sue  rending  on  a  Quartz 
Spring  in  an  Evacuable  Chamber 
This  method  is  a  traditional  technique  for  sorption 
measurements.  Figure  (20)  illustrates  the  essential 
components  for  this  system.  The  chamber  is  evacuable  and 
under  thermal  control.  It  connects  with  a  vapor  source  to 
cause  a  mass  change  of  sample  which  can  be  measured  through 
a  telescope.  Figure  (21)  shows  the  main  parts  in  more 
graphic  detail.  Figure  (22)  depicts  the  vacuum  system.  The 
thermal  scarce  contains  a  50  Q  heater,  a  temperature 
controller  model  49  (or  50),  a  10  switch  OMEGA's  digital 
temperature  indicator  (model  175),  a  centrifugal  fan,  and  a 
variator.  Six  thermocouples  attached  no  the  system  allow  us 
to  calibrate  the  temperature  of  different  parts  of  the 
system.  The  curves  shown  on  Figure  (23)  and  (24)  show  these 
temperatures  and  will  be  discussed  below.  From  Figure  (24), 
it  is  easy  to  see  the  time  needed  to  reach  steady-state  is 
a:  least  50  rain.  The  vapor  source  is  a  removable  tube 


components  in  the  apparatus  for 
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Fig. 22  Vacuum  system 


TEMPERATURE  CALIBRATION 


Fig.  23 
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connected  to  the  chamber  through  a  three-way  stopcock.  The 
vacuum  system  is  composed  of  a  combination  of  a  mercury 
diffusion  pump  and  a  mechanical  pump. 

In  order  to  demonstrate  the  reliability  of  this  set 
up  PET  (Poly-ethylene  Teraphtnalate)  film  was  chosen  because 
extensive  reliable  data  for  H20  diffusion  in  this  polymer 
were  available.  The  sample  was  10.15  cm  long,  1.57  cm  wide 
and  80  um  thick.  Obviously  the  edge  effect  is  insignificant. 

(See  Appendix  D) .  Figure  (25)  gives  sorption  and  desorption 
curves  of  PET  at  36°C.  Figure  (25)  gives  the  curves  for 
sorption  at  three  different  temperatures.  Diffusion  co¬ 
efficients  were  calculated  based  on  the  initial  gradients  of 
the  curves.  The  slopes  were  auite  linear  up  to  M^M  ~  60*. 

*  u  ^ 

This  is  called  Fickian  behavior.  Figure  (27)  is  an  Arrhenius 
ploe  of  H^O  diffusion  in  PET.  The  result  shows  good  agree¬ 
ment  with  literature  values. 

The  activation  energy  10.4  kcal/mol  (43.5  kJ/mol)  is 
the  same  as  is  given  in  the  literature.  The  deviation  of 
diffusion  coefficients  are  within  10-5.  These  results  demon¬ 
strate  that  the  system  is  reliable.  That  is  not  to  say  than 
problems  were  not  still  there.  For  example  the  tiny  weight 
of  the  PP3T  samples  was  beyond  the  sensitv  of  our  original 
quarts  sprma ,  then  a  new  quartz  spring  with  10  mg  maximum 
load  and  47.5  ag/mm  sensitivity  came  to  replace  the  old  one. 


AT  36°  C 


i> 
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which  was  200  mg  maximum  load  and  9.75  mg/cm  sensitivity. 

The  new  spring  is  so  fine  that  it  is  almost  impossible  to 
handle.  To  my  knowledge,  the  manufacturing  company  can  not 
make  any  finer  quartz  springs  than  this. 

One  protecting  column  was  added  to  the  system  to 
protect  the  spring  from  being  hit  by  air  or  vapor  flow,  or 
being  sucked  into  the  vacuum  system.  The  solubility  of 
diffusar.t  in  PPBT  came  out  to  be  so  low  than  the  thermal 
expansion  of  the  chamber  had  to  be  considered  in  some 
experiments.  The  sample  used  in  this  method  is  23555-25-6 
PPBT  film  with  weight  9.658  mg,  length  ~5.2  cm, ' width  ~6  mm, 
and  thickness  20  pm. 


Chapter  IV  RESULTS  AND  DISCUSSION 


t  of  Microcavities  on  Diffusion 


I  have  previously  mentioned  that  the  PPST  film  avail¬ 
able  for  these  experiments  is  still  very  rough.  The 
surface  contained  erupted  voids  which  can  be  seen  at  high 

magnification  in  the  electron  microscope.  According  to  a 

37 

report  from  the  University  of  Massachusetts,  they  rer erred 
to  the  microvoids  inside  PPBT  fibers  and  films  as  needle¬ 
like  in  shape.  Perhaps  football-like  would  be  an  alterna¬ 
tive  description.  The  Average  Void  Sizes  obtained  by  Debye 

O 

analysis  of  small  angle  X-ray  scattering  is  163  ±  10  A 

O 

along  the  extrusion  direction  and  135  i  11  A  perpendicular 
to  the  extrusion  direction.  This  correlates  to  ~13  monomers 

O 

(12.35  A  each)  along  the  extrusion  direction  and  ~11 

O 

monomers  (still  based  on  length  12.35  A)  perpendicular  to 
the  extrusion  direction.  These  microvoids  will  act  as  deep 
trapping  sites  in  diffusion.  When  penetrants  are  sorbed 
into  voids,  they  become  virtually  immobilized.  Naturally, 
the  observed  diffusion  rate  will  decrease. 

Figure  (28)  shows  three  curves  for  sorption  rates 
measured  by  the  electrobalance.  The  changing  of  room  con¬ 
ditions  and  the  different  time  delays  for  transferring 
samples  from  the  desiccator  to  the  electrobalance  correspond 


co  the  main  reason  for  the  curve  shift. 


Another  factor  is 


the  drifting  of  the  electrcbaiance ,  which  has  been  mentioned 

-10  "> 

before.  After  calibration,  the  D  value  is  6.0  :<  10  cm“/ 

-3 

sec  (see  Figure  (29))  and  solubility  is  1.28  x  10 

Compare  these  values  to  the  data  from  the  quartz  spring 

-10  2  - 
mechod  which  has  D  =  1.79  x  10  cm  /sec  and  s  =  5.23  x  10 

(See  Figure  (30),  (31))  The  results  are  dramatically 

different.  One  is  forced  to  inquire  "What  causes  this 

difference?"  There  are  some  important  differences  in  the 

way  the  experiments  were  done.  The  experiment  v/ith  the 

electrobalance  occurred  in  air  with  70%  humidity.  On  the 

other  hand,  the  experiment  in  which  the  quartz  spring  was 

used  is  in  vacuum  initially  and  then  has  almost  100% 

humidity.  But  these  factors  would  not  be  expected  to  make 

any  big  difference  in  the  measured  diffusion  rates  for 

normal  rubbery  polymers  or  other  homogeneous  polymer  systems 

Due  to  the  fact  that  PPBT  has  lots  of  microcavities,  a 

26 

"dual  sorption  theory"  was  introduced  which  enabled  us  to 
resolve  the  conflict.  First,  we  assume,  when  the  sample 
is  measured  in  air,  that  the  microvoids  have  been  filled 
with  air.  Water  will  not  be  trapped  in  these  voids.  Then, 


we  apply  equation  (66)  from  chapter  II, 
D_^  =  D/1  +  R 


ycr.i  =  1.79  x  10  cm-/ sec  that  is  the  value  frorr.  the 

quartz  spring  and  D  =  6.0  x  10  crr.^/sec  which  is  the 
value  from  electrcbalance . 

We  can  get 


=  2.35,  where  R  = 


bS,'  /K , 
h  a 


bS'  =  2.35K, 
h  a 


In  this  case,  the  pressure  is  low  sc  equation  (53)  is  used 


C  =  (Kd  +  bS^Jp 


=  (K,  +  2 . 3 5 )  p 


=  3 . 3  5K  ,p 


where  C  =  5.23  x  10  which  is  the  solubility  of  M 0 0  in 
using  the  method  of  quartz  spring.  Therefore, 


K^p  =  1.56  x  10 


According  to  our  assumption  that  water  won't  be  trapped 
inside  the  microvoids  by  using  the  electrobalance.  The 
solubility  will  follow  Henry's  law.  So,  the  value  can  be 
calculated  by  putting  Q.7p  (because  of  7GG  R.H.  in  room) 
into  above  equation. 


Kdx(0.7p)  =  1.09  x  10 


-3 


Now,  comparing  this  value  to  the  experimental  result, 

-3 

1.23  x  10  from  the  electrobalance,  the  results  are  in 

reasonable  agreement  with  each  other.  A  little  bit  high' 

value  from  electrobalance  implies  that  the  microcavities 

still  trapped  sene  small  quantity  of  water. 

From  the  above  demonstration,  the  water  sorbed  by 

microvoids  at  25 5C  is  approximately  3.57  x  10  That 

corresponds  to  3.54  x  10  “  mg.  If  the  water  we re  in  the 

-5  3 

liquid  phase,  the  volume  should  be  3.54  x  10  cm  . 
Dividing  by  the  sample  volume  {6.24  x  10-J  cmJ)  yields  ti 
volume  ratio  5.67  x  10  ^ .  Obviously,  the  water  is  not 
saturated  in  the  voids  unless  the  vapor  pressure  is  very 
high.  Consequently,  the  voids  volume  should  be  much  big: 
than  this  amount.  That  is  voids  volume  5.57  x  10  ^ 
(0.567  v/o) . 

3.  Activation  Energy  Relates  to  Rotation  Energy  of  Chai; 
Segment 


P?3T  film  is  anisotropic  in 
Hsu  et  al.  of  the  University  of 
wide  angle  electron  diffraction 
prooable  monoclinic  unit  cells 
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c 1  =12. 35  A  c  =  12.25  A 

Y  =  96°  y  =  63° 

2  =  1  z  =  2 

The  proposed  arrangement  of  PPBT  in  its  crystal  structure  i 

shown  as  Figure  (32) .  How  does  this  relate  to  diffusion? 

According  to  the  molecular  model  proposed  by  DiBanedetto 
°0 

and  Paul  the  activation  energy  for  diffusion  is  related 
to  the  variation  in  the  average  potential  energy  of 
molecular  interaction  which  accompanies  a  volume  change  of 
\n(n/4)r5\ 
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Pace  ar.d  Dauyner,  following  this  idea,  proposed  a 
statistical  mechanical  model  in  which  the  activation  energy 
is  related  to  the  bending  energy  and  Van  der  Waals  energy 
for  producing  a  separation  d  between  the  surfaces  of 
adjacent  chains  where  d  is  the  diameter  of  r.he  diffusant. 
Both  theories  point  out  that  the  activation  energy  for 
diffusion  is  the  energy  necessary  to  produce  and  effective 
volume  large  enough  to  accommodate  the  diffusant.  Now,  we 
are  dealina  with  a  polymer  whj.cn  has  flat  rigid  chain 
segments.  (See  Figure  (33)).  It  displays  the  possibility 
of  rotation  on  the  two  single  bonds. 


According  to  the 


Cl 


segment  ~o 


riecnam; 


roxatiai;  me  cr.a;n 


aecoaoca 


accomodate  benzene  molecule. 


report  or  J.  E.  Mark,  et  al.  of  the  University  of 
Cincinr.atti,  that  the  ring  groups  freely  rotate  to  20°  from, 
the  origin.  The  activation  energy  is  2  kcal/mol.  If  it 
rotates  to  90°,  the  rotation  energy  is  6  kcal/mol.  The 
activation  energy  for  I^O  diffusion  in  PPBT  is  5  kcal/mol 
and  that  for  benzene  diffusion  is  12  kcal/mol.  (See  Figure 
(31),  (34)).  This  implies  a  diffusion  mechanism  as  suggested 

in  Figure  (33).  H^O  needs  only  one  chain  segment  to  rotate 
90°  for  it  to  pass  through;  however,  benzene  needs  both 
chain  segments  to  rotate  90°  for  it  to  get  through.  If  one 
remembers  the  discussion  in  the  previous  section  of  how 
microcavities  affect  diffusion,  it  can  be  seen  that  the 
activation  energies  measured  are  not  the  true  values.  The 
effective  diffusivity  Dgf -  =  D/1  +  R  will  change  with  R. 

If  R  decreases  when  temperature  increases,  the  measured 
effective  activation  energy  will  be  higher  than  it 
supposed  to  be.  Generally,  this  is  the  tendency  that 
dlnb/d(l/T)  higher  than  dlnK^/d (i/T) .  Consequently,  this 
will  cause  an  abnormal  high  negative  enthalpy  of  solution. 
Figure  (35)  and  (36)  show  the  evidence  for  this. 


C.  Interpretation  of  Coefficients 

The  sample  of  PPBT  has  a  very  rough  surface,  uneven 
thickness,  voids,  and  impurities.  Rough  surface  and  uneven 


Kv'CviSis 


IENTUS 


000/T 


SOLUBILITY  COEFFICIENTS  OF 
WfiTER  IN  PPBT 


thickness  can  even  be  seen  by  the  naked  eyes.  Voids  nave 
been  studied  by  using  snail  angle  X-rav  scattering  patterns. 
That  voids  exist  in  substantial  amounts  also  has  been 
supported  by  previous  discussion  in  this  thesis.  The 
probable  impurities  in  PPBT  are  residual  methane  sulfonic 
acid  (or  PPA  which  has  been  found  in  2DAX  in  this  thesis) 
and  other  impurities  which  could  be  introduced  in  the 
synthesis  procedure.  All  these  effects  will  shift  the 
diffusion  character  significantly  from  the  theoretical  pre¬ 
dictions  of  well  dimensioned  pure  samples.  Table  I  gives 
the  coefficients  of  water,  benzene,  and  ethanol.  The  lack 
of  enough  experimental  data  on  ethanol  is  because  the 
solubility  of  ethanol  is  too  low  to  be  detected  at  high 
temperature.  The  temperature-dependence  of  the  solubility 
coefficient  includes  the  heat  of  condensation  and  the  heat 
of  mixing.  In  our  case,  the  Lane  iuir  isotherm  of  solution 
in  pre-existing  cavities  also  exists.  The  abnormal  enthalpy 
of  solvation  is  mainly  attributed  to  the  S^b  quantity  which 
is  the  hole  filling  term  in  S  =  (K^  +  bS^)p  as  discussed  in 
chapter  II.  The  order  of  solubility  coefficients  in  terms 
of  mole  numbers  is  water  >  benzene  >  ethanol.  In  contrast, 
the  diffusion  coefficients  at  room  temperature  have  the 
order  as  ethanol  >  benzene  >  water.  If  we  recall  that 
D  f  -  =  D/(l  +  (bSf  /  K , )  )  and  S  =  (K,  +  bS,')p,  it  is  obvious 


Water 


Ethanol 


Ed (kCal/mole) 


2 

D  (cm  /sec) 


H  (kCal/raola) 
s 


2 

DQ(cm  /sec) 


S (mg/mg) 


Ben cone 


6 


12 


_i  o  -io 

1.71x10  2.8x10 

(25  °C)  ( 24  °C) 


1.1x10 
( 24  °C ) 


-9  (high  T)  -23  ^ 

-32  (low  T) 


4.28x10 


0.188 


a. 


6.93x10 
( 25  °C ) 

(24+1 . 4mmHg) 


7.99x10  3  1.48x10  3 

( 24  °C )  ( 24  °C ) 

(94+4mmHg)  (55±lminHg) 


Table  I 

Diffusion  Data  of  PPBT  (28555-25-6)  Film. 


that  the  higher  the  bS^  value  the  higher  the  S  value,  but 

the  lower  D  value.  If  we  further  assume  that  SJ  is  the 
ef  f  h 

same  for  these  three  species,  it  will  imply  that  the  forces 
between  polymer  surface  and  diffusants  are  in  the  order 
water  >  benzene  >  ethanol.  Figures  (37)- (43)  are  diffusion 
curves  for  the  three  vapor  species  studied.  The  solid 
curves  represent  theoretical  values .  All  data  display  a 
negative  deviation  in  the  long  time  region.  This  is  the 
behavior  of  a  concentration  dependent  diffusion  rate.  Again 
microcavities  are  responsible  as  will  be  demonstrated  below. 
The  flux  of  penetrant  is  only  due  to  mobile  molecules,  not 
trapped  ones,  therefore 


J  =  -D  (C) 3C/3x 


-  -D1 (C) 3C1/; 


where  C  is  the  overall  concentration,  the  suffix  1  refers 
to  mobile  molecules.  We  find  D(C)  =  (C) 3C^/ 3C .  If  the 

immobilized  molecules  are  sorbed  according  to  a  Langmuir 
model  and  the  mobile  molecules  according  to  Henry's  law, 
then 


C2  =  Cl/  (ct  +  0C1) 


vnore  C„  denotes  the  concentration  of  immobile  molecules, 


ifi.  38 


3C/5C,  =  a/  (a  +  3th) 


d  r*  C*  30 


d.  j;  -  D^Oj/11  -  ec) 


This  simple  model  leads  tc  a  diffusion  coefficient  which 


increases  with  concentration. 


Thus  ,  the 


:ed  materials 


acting  as  traps  inside  PPBT  gi\re  the  same  effect.  This 
phenomenon  also  was  observed  in  this  experiment. 

Figure  (44)  displays  the  sorcerer,  ana  desorption 
curves  of  H-,0  in  PPBT.  The  curves  are  sigmoid  in  shape 
while  the  desorption  curve  crosses  the  sorption  curve  twice 
before  equilibrium  is  reached.  This  sort  of  behavior  which 
is  described  as  anomalous  or  "non-Fickian "  has  been  re- 


portec , 


A.  —  —  .4.  ±  A  _> 


:ar.ce,  in  the  sorption  of  rater  by  cellulose 


39  4  0 

(Mewns,  1956),  Keratin  (King,  1945)  and  vinyl  alcohol 

4  1 

(Long  and  Thompson,  1955). 

The  water  remaining  inside  PPBT  in  the  desorption 
procedure  is  1.02  :<  10  ^  mg/mg. 


As  mentioned  in  chapter  II  (part  G)  relative  to  ionic 


conductivitv 


water  will  enhance  the  dissociation  of  ionic 


mpuricios  into  anions  ana  cations,  and  consequently  m 
he  conductivity.  Will  any  impurities  in  polymers  aff a 
if fusion?  Yes,  this  is  demonstrated  by  measurements  o 

ms  soaked  for  2  days 
1  molar  salt  solutions  of  three  nitrates  (Fe(NO. )  -  , 

j  a 

Co  (NO^ )  2  /  NifNO^)-,)  and  two  chlorides  (CaC^,  LiCl)  .  PPB 
samples  were  immersed  in  1  molar  salt  solutions  for  two 
days  and  rinsed  in  distilled  water  before  being  dried.  T 
dried  samples  were  measured  by  the  electrobalance .  From 
the  rates  at  which  the  weight  increased,  diffusion  co¬ 
efficients  were  determined  (Figure  (45),  (46)).  A  plot  o 

logD  vs  the  charge  to  radius  ration  q/r  came  out  tc  be  a 
straight  line  (Figure  (47)).  This  may  be  described  by  th 
following  equation 


logD  -  0.214q/r  -  3.592 


where  q  =  the  electric  charge  of  cation 
r  =  crystal  ionic  radius 

The  constant  -9.592  corresponds  to  D  =  2.56x10  ^ 
cm“/sec.  This  value  is  very  close  to  the  diffusion  co¬ 
efficient  of  water  in  pure  PPBT  using  the  quartz  spring, 
suggests  that  the  ions  are  mainly  in  microcavities. 

Lot  us  consider  any  rate  constant  K  can  be  expressed 

—  F  /  RT 

as  an  Arrhenius  equation  K  =  Ae  ,  then  InK  =  InA  -  E/ 


•An 


comcare  this  to  equation  (78),  we  find  that  if  water  acts 


on  another  neaative  charae  -q  then  -qc  /r  correlates  to 

-  w  *w 

the  energy  term  E .  Thus,  the  Diffusion  rate  is  attributed 

to  the  reaction  between  water  and  cations.  Also,  the 

diffusion  rate  is  directly  proportioned  to  the  electrical 

potential  fields  of  cations.  One  experiment  cone  by 
49 

Takamatsu  about  sorption  phenomena  in  nation  Membranes 
is  described;  during  the  neutralization  process  which 
cations  replaced  H+  in  SO^  H+  site,  the  coefficients  D  of 
several  cations  were  determined  and  logD  was  related 
linearly  to  q/a,  where  a  is  the  separation  distance  between 
centers  of  charge  of  the  cation  and  anion.  This  observa¬ 
tion  might  have  the  same  reason  as  ours. 


E.  Possible  Interaction  of  Doped  Salt  with  Polymer 

There  is  a  very  important  matter  for  organometallic 

A 

chemistry,  called  the  "sixteen  and  eighteen  electron-rule." 
According  to  this  rule,  stable  organometallic  compounds  of 
the  transition  metals  will  have  a  total  18  valence  electrons 
about  the  metal;  in  other  words,  they  will  have  the 
"effective  atomic  number"  (EAN)  of  the  next  higher  inert 
gas.  Molecules  having  only  16  valence  electrons  can  often 
be  just  ns  stable  as  13".  In  general,  organic  complexes  of 
the  transition  metals  are  formed  with  those  metals  in  low 


v,  >V> 


11 


oxidation  states,  since  it  is  only  in  that  condition  chat 

the  metal  will  have  populated  orbitals  of  a  symmetry  and  a 

* 

low  z  to  allow  metal  to  ligand  a  electronic  flew.  Thus, 

2+  •+■  0 

the  prec erred  oxication  stares  are  Fe  ,  Co  ,  Ni  .  Un¬ 
fortunately,  these  are  not  the  ones  in  our  case.  But, 
possible  reactions  will  happen,  if  reduction  of  metal  ions 
occurs . 

Diffusion  rates  were  measured  repeatedly.  (Figure 
(43) -(30),  Tables  II-V) .  Since  the  room  conditions  changed 
and  due  to  the  limited  sensitivity  of  instrument,  it  is 
still  not  possible  to  make  final  conclusions.  However,  ic 
is  obvious  that  the  diffusion  rates  decreased  to  some  exte.n 
The  sample  weight  decreased  somewhat  ana  in  some  cases, 
even  appeared  to  become  lower  than  that  of  the  original  dry 
sample • 

F.  Discussion  on  .’orpnology  and  Anisotrouic  Charac car . 

If  diffusion  in  PPBT  is  a  process  of  mass  transport 
along  the  interchain  channels,  then  the  diffusion  rate 
should  be  sensitive  to  the  morphology  of  the  polymer.  It 
is  then  necessary  to  analyze  the  relation  between  morcnolog 
and  diffusion. 

2 

According  to  Edwin  L.  Thomas  et  al.,  the  PPBT  tiber 
can  be  peeled  and  mildly  sonicated  into  thinner  and  thinner 
ribbonlike  fibrils.  Each  of  these  fibrils  appears  to 
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Taole  II 


Sorption  Coefficients  of  H2O  in  Nitrate  Doped 
PPDT  Pi  Ini  After  Several  Days  (before  correction). 
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Table  TV 

O 

Diffusion  Coefficients  (cn'Vs)  )i  Water  in  Nitrat 
Doped  PPBT  Film  After  Several  Days. 


Pure  dry  Samp  la 

befora  doped  6.36  7.11 

(mg) 


1st  day 

doted  6.901  7.16 

(mg) 


2nd  day 

doped  6.38  7.16 

(mg) 


3rd  day 

doped  a,  7.16 

(mg) 


4th  day 

doped  6.86 

(mg) 


30th  day 

doped  6.82  7.11 

(mg) 


Table  V 


Weight  (mg)  of  PP3T  Film  Measured  Before 
After  Doping  for  Several  Days  of  Drying. 


cc r.sist  or  sn’.ri.er  micro: icr ils  or  lateral  dimensions 

o 

varying  from  50  to  30  A.  (Figure  (51))  .  A  schematic  of  t 
texture  of  these  kinked  ribbons  is  proposed  in  Figure  (52) 
Each  microfibril  consists  of  a  succession  of  narrow 

O 

crystallites  (50  x  500  A  or  smaller)  embedded  in  somewhat 
less  ordered  matrix.  The  periodicity  of  the  banding  varie 

O 

from  1000  to  2000  A.  At  present,  they  still  do  not  know 
whether  the  bands  that  appear  are  due  to  the  fragmentation 
of  the  fibers  during  the  sample  preparation,  or  on  the 
contrary,  are  characteristic  of  the  as-spun  fibers. 

Based  on  above.  Barker  proposed  that  microvoids  might 


exist  between  microfibrils  as  suggested  by  Figure  (53. a.). 
Then  the  parallel  diffusion  will  be  much  larger  than 

O 

transverse  diffusion.  If  the  voids  are  only  a  few  A  width 
the  longitudinal  diffusion,  parallel  to  the  axis  of  micro¬ 


fibrils,  and  the  activation  energy  for  a  jump  x  is  approxi 
matelv  proportional  to  a  quantum  y  of  the  order. 


*  *  y 

:  P  /a 


From  equation  Df/  =  gL“  ■;  ,  wc  have,  as  a  first 
aoproximation 


x/-d:h/ffe--:/RTcix) 


2  2 
!RT)  ! 


Insert  some  assumed  values  to  above  equation,  (See  appendix 
E) 

a  -  1.235  A,  e*  ~  1645  Jmol-1 

*  °  12-1 
p  ~  4.3  A,  Vq  =  10  sec 

-5  2 

it  gives  the  result  Dy/  ~  2  x  10  cm  /s  or  the  order 
-10  ^  cm“/s  at  temperature  300°K.  If  the  microvoid  is  much 
bigger  than  a  few  angstroms;  diffusion  rate  will  approach 
to  gas  diffusion  in  air,  for  example,  water  vapor  is 
0  -  0.239  cm^/s  at  3°C. 

The  microvoids  might  not  exist  if,  instead  of  voids, 
the  intervals  were  filled  with  less  ordered  matrix.  In 
this  case,  Dy/  could  be  only  a  few  times  greater  than  D± . 
Actually  the  present  available  PP3T  samples  were  demon¬ 
strated  to  have  lots  of  microvoids  of  large  sizes: 

O  O 

(163  t  10  A  135  ±  11  A  for  heat  treated  film  and  49.1  ± 

O  O 

3.5  A  x  34.7  ±  2.5  A  for  heat  treated  fiber)  as  mentioned 
in  section  one  in  this  chapter.  Despite  this  fact,  it  is 
still  interesting  to  extend  the  theoretical  approach. 

To  simplify  the  question,  we  assume  the  PPBT  is  1003 
crystalline.  Then  diffusion  occurs  only  when  the  diffusant 
can  penetrate  into  the  crystal.  Under  this  assumption,  the 
approach  in  section  2.  of  this  chapter  can  be  applied  to 
the  calculation  of  the  ratio  Df/  /D^  .  The  proposed  parallel 
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diffusion  mechanism  is  displayed  in  Figure  (53.b.).  Recall 
from  equations  (44)  and  (45),  the  A  value  should  be  the 
same  for  both  cases.  Thus 

D7/  /D_  =  L^/LL2exp(- (AE/,  -  AEJ  /RT) 

—  O  _ 

is  between  3.5  and  6.1  A  and  L/,  is  half  of  the  value 

O 

12.35  A,  because  the  rotating  segment  will  be  either  the 
phenyl  group  or  the  bencobisthiazole  group. 

—  O  _  O 

Assume  L,  is  3.5  A,  L//  is  6.13  A, 
if  AE/(f  -  AE±  =  -4  kcal/mol, 
then  D/y  /D^  ~  2.67  x  IQ2. 

If  AE y.  -  AEX  =  0  kcal/mol, 
then  D t/  / D.  ~  3.1. 

From  above,  it  is  easy  to  see  that  D/(,  /D,  is  mainly 

dependent  on  the  difference  between  AE ,  and  AE , . 

3efore  ending  this  chapter,  I  ’would  like  to  point  out 

that  the  viscosity-diameter  of  a  diffusant  is  usually 

0  21 

large.  For  example,  0^  has  a  d-value  3.5  A.  Therefore, 
it  is  impossible  to  accommodate  a  diffusant  without  further 
bending  after  rotation.  The  discussion  in  section  2.  of 
this  chapter  was  based  on  molecular  shape.  Also,  the 
protor.ated  P?PT  will  have  some  tilt  angle  on  its  chain 
segments,  which  is  not  shewn  on  our  Figure  (53.b.)  and  our 
discussion  has  not  dealt  with  that  for  the  purpose  of 
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Chapter  V  CONCLUSION 

It  is  proposed  that  the  existence  of  microcavities  in 
PPBT ,  lead  to  the  measured  diffusion  rates  which  were 
smaller  than  the  true  values.  The  evidence  for  the  micro¬ 
void  model  are:  1.  There  is  an  abnormally  high  heat  of 
sorption  (solution).  2.  There  is  a  lower  water  solubility 
when  exposed  to  water  in  air  than  when  exposed  after  being 
in  a  vacuum.  3.  The  diffusion  rates  are  concentration 
dependent.  These  phenomenon  fully  consistent  with  the  dual 
sorption  theory. 

The  expected  high  activation  energies  were  not  found. 
This  is  thought  to  be  due  to  the  rotation  or  partial  rota¬ 
tion  of  chain  segments  on  the  mer.  Statistical  mechanical 
model  of  Pace  etc.  is  a  powerful  tool  for  predicting  activa¬ 
tion  energy  with  its  direct  calculation  using  non-ad justable 
parameters.  It  would  be  worthwhile  to  make  detailed 
theoretical  calculations  if  enough  precise  data  were  avail¬ 
able  . 

The  sigmoid  shape  of  sorption  and  desorption  curves 
and  negative  value  of  logD^  in  the  case  of  water  diffusion 
are  typical  of  polymers  below  their  glass  transition  tempera¬ 
ture.  Actually  the  PPBT  is  in  crystalline  state. 

The  diffusion  rates  of  H^O  in  salt  doped  PPBT  are 
linearly  proportional  to  q/r.  This  effect  is  consistent 


with  the  idea  that  the  cations  act  as  traping  sites  by- 


attracting  water  with  an  electrical  field. 

The  evidence  that  doping  salt  is  deposited  mainly  in 
microvoids  are:  1.  D  at  q/r  =  0  has  the  value  close  to 
that  of  pure  PP3T  measured  subsequently  to  being  in  the 
vacuum  chamber.  2.  The  sorption  coefficients  cf  H^O  in 
Co(hT0-J9  and  in  NitNO^)^  doped  PPBT  are  approximately  equal 
to  the  value  measured  of  pure  PPBT  in  the  vacuum 
chamber  method.  Interactions  between  doping  salt  and  PPET 
are  not  conclusive  but  significant  changes  in  the  diffusion 
curves  and  solubilities  are  obvious  after  a  long  time  with 
the  salt  ir.  polymer. 


r. 
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APPENDIX  A. 


ERROR  FUNCTION  SOLUTION  TO 


DIFFUSION  EQUATION 


The  equation  of  linear  flow  mass  is 


D82C/3x“  =  3C/3t 


(79 


Consider  the  expression 


r  _  ,  -^^-xi'/4Dt 
C  —  t  c 


Since 


3C/3t  =  -  (l/2t3/2)e~X  /4Dt+  (x2/4Dt5/2)e  x  /4D< 


and 


2  2 
3  C/3x 


-  (l/2Dt3/2)e_x  /4Dt-r  (x2/4D2to/2)e  X  /4Dt 


It  follows  that  equation  (SO)  is  a  particular  solution  of 
euqation  ( 79 ) • 

Clearly,  further  solutions  of  equation  (79)  are 
obtained  by  dif ferentialing ,  or  in  some  cases  by  inzergrat 
ing,  equation  (80)  with  respect  to  either  x  or  t.  So,  tha 
equation  (79)  is  satisfied  by 


x  2,,  u  x/2(Dt)'4  _ r 2 

C  4-  J  X  /lUt,  _  f  , , 


The  notation 


X  2 

erf  x  =  (2//tt)J  e  ^  d£ 

0 

will  always  be  used  here,  so  we  have  just  shown  that 
A  erf  (x/2  /Dt  ) 

whore  A  is  an  arbitrary  constant,  is  a  solution  of 
equation  (19) . 

7  2 

Because  D5“C/3x  =  3C/3t  is  a  second  order  partial 

differential  equation  in  x,  the  solution  can  be  expressed 

C(x,t)  =  A  erf(x/2/Dt)  +  lx  +  m  (3 

New,  the  boundary  conditions  are  introduced  to  the 
equation  as  below: 


C(x,t)  =  C^,  x  =  0  for  all  t, 

C(x,t)  =  Ch,  x  >  0,  t  =  0. 


According  to  the  boundary  conditions  above,  when  t  =  0, 
x=  0,  C(x,t)  is  a  constant.  So,  the  constant  1  should  be 


equal  to  zero,  then  equation  (31)  becomes 
c(x,t)  -  A  erf  (::/2,'CT:)  f  m 


Co  cause 


erf(O)  =  0,  erf  (■»)  =  1 

by  applying  the  boundary  conditions  to  equation  (32) 
have 

m  =  C 

s 

A  +  m  =  C . 

l 

Consequently 

A  =  C.  -  C 
i  s 

Therefore 

C(x,t)  =  (C.  -  C  )  erf  (x/2/Dt)  +  C 
is  s 

+  C  (1  -  erf  (x/ 2 * ITt) 

+  C  eric  (:</  2  v  Dt ) 


=  erf (x/ 2 /Dt) 
=  C.  erf (x/2/Dt) 


APPENDIX  B  THE  SOLUTION  TO  DIFFUSION 


EQUATION  WITH  INHOMOGENEOUS 


BOUNDARY  CONDITIONS 


The  boundary  conditions  are  as  below: 


C(x,t)  =  at  x  =  0,  all  t  >  0 


C(x,t)  =  C2  at  x  =  b,  all  t  >  0 


C(x,t)  =  Cq  at  0  <  x  <  b,  at  t  =  0 


where  b  is  the  thickness  of  polymer  membrane,  and  C0 
are  the  concentrations  of  gases  just  inside  the  polymer 
surfaces  and  Cq  is  the  initial  concentration  of  diffusant, 
assumed  uniformly  distributed. 

In  order  to  solve  the  equation,  one  new  concentration 
variable  U(x,t)  is  introduced,  i.e. 


u(x,t)  =  c  (x ,  t )  +  n(x) 


C(x)  is  a  first  order  linear  function  of  x,  because  the 
2  2 

equation  3  ft(x)/3x  =  0  must  be  satisfied  for  our  purpose  of 
solving  the  equation.  It  is  defined  as  C(x)  =  ax  +  8 • 

Now,  the  new  boundary  conditions  are  set  as: 


U (x, t)  =  0  , 

at  x  =  3  for 

all 

t 

U(x,t)  =  0, 

at  x  =  b  for 

t  = 

0 

(84 

U(x,t)  =  C 

+  ft  (x) ,  at  0 

<  X 

<  b,  at  t  =  0 

Pick's  second  law  becomes: 

3U(x,t)/?t  =  D(  U2U(;:,b)/jX2)  +  ( 32ft (x) / 3x2) } 

=  D  (92U(x, t) / 3x2 } 

Here  U(x,t)  =  C(x,t)  +  ft(x)  must  satisfy  the  boundary 
condition  (34) .  That  is 

at  x  =  0,  +  3  =  0,  then  3  =  -C^ 

at  x  =  b,  C9  +  ab  +  8  =  0,  then  a  =  (C^  -  C2)/b 

Therefore , 

U(x,t)  =  C  (x ,  t )  +  (  (Cx  -  C2)/b)x  -  Cx 


(x,  c)  -  U 


vx  , 


-  (  « 


:.)/b)x 


(85 


We  begin  our  attack  on  the  problem  by  attempting  to 
lir.a  solutions  of  the  diffusion  equation  that  are  of  the 


soncial  form 


U  ( x  ,  t )  =  X  (x)  T  (t ) 


(3 


We  impose  the  requirements 


X(0)  =  0,  X  (b)  =  0 


in  order  that  these  "product  solutions"  will  also  satisfy 
the  boundary  conditions 


U(0,t)  =  U  (b ,  t)  =  0 


If  the  diffusion  equation  has  a  solution  of  the  form 


(36),  then  we  must  have 


X  (x )  T  '  (t)  =  DX”  (x)T(t) 


or  upon  dividing  through  by  DX(x)T(t) 


X "  ( x )  / X  ( x )  =  -T'(t)/DT(t) 


The  left  hand  side  of  this  equation  is  independent  of  t 
while  the  right  is  independent  of  x.  Therefore,  both 
members  must  be  equal  to  a  constant,  which  we  denote  by 
thus 


X "  ( x ) / X ( x )  =  -T'(t)/DT(t)  =  -A 


and  so  the  function  X  and  T  must  be  solutions  of  the 
ordinary  differential  equations 


-i.  L  .<W 


mvmmw 


are  solutions  of  th 


respectively.  Conversely,  if  X  and  T 
corresponding  equations  for  the  same  value  of  A,  then  the 
product  U(x,t)  =  X(x)T(t)  is  a  solution  of  the  diffusion 
equation  because 


U.  (x,t)  - 


DU  ( x , t ) 

XX  ' 


=  X  { X ) T 1  (t)  -  DX" ( x ) T { t ) 


=  X(x)  (-\DT(t)  )  -  D(-,\X(x)  )T(t) 


=  0 


The  requirements  for  X,  namely 


X"  +  XX  =  0,  X (0)  =  0,  X (b)  =  0 


constitute  an  eigen  value  problem.  The  eigen  values  are 


Xn  =  (n-rr/b) 


n  =  1,2,3,... 


The  corresponding  eigenfunctions  (non  trivial  solutions) 
are 


Xn  =  sin  (nrrx/b) 


n  —  j.  ,  2 , 3  ,  .  .  . 


■or  X  =  \  ,  we  have  T '  +  }  t  =  0 .  A  solution  i: 
n  n 


Tn(t)  =  exp ( - X^Dt ) 


=  exp(-n“a  Dt/b“) 


Each  of  the  product  functions 


Un(x,t)  =  exp  (-n2n2Dt/h‘')  sin  (nnx/b)  ,  n  =  1,2,3,... 

is  a  solution  of  the  diffusion  equation  and  satisfies  the 

boundary  conditions  (87) .  If  the  constant  C  are  such  that 

n 

the  series 

CO 

u(x,t)  =  l  C  U  (X,t) 
n=l 

CO 

=  7  C  exp  (-n^T^Dt/b^)  sin  (n;rx/b)  (38) 

n=l  n 

converges  and  can  be  differentiated  termwise  a  sufficient 

number  of  times  with  respect  to  x  and  t,  then  this  series 

also  represents  a  solution  of  the  differential  equation 

that  satisfies  the  homogeneous  boundary  conditions  (87) . 

However,  the  constants  C  must  be  chosen,  if  possible, 

n 

in  such  a  way  that  the  series  (88)  satisfies  the  non- 
homogeneous  boundary  condition 

U  ( x  ,  0 )  =  f  ( x )  ,  0  <  x  <  b 

Thus,  we  require  that 

CO 

U(x,0)  =  y  C  sin  (mx/b) 

J  i  n 

n=l 

=  f(x)  0  <  x  <  b 


s 


"i 


But  tha  C  must  be  the  nth  coe::_cien:  m  tha  Founsrsir.a 
n 

series  for  f,  that  is 
b 

C  =  (2/b)  /f  (x)  sin  (mrx/b)  dx  (39) 

n  A 


Here,  f (x)  =  CQ  +  fl(x) 

=  C0  -  C1  +  ( (C1  -  c2) /b) x 

can  be  used  to  find  C  . 

n 

From  equation  (89)  ,  we  have 
b 

Cn  =  (2/b) J  (CQ  -  C1  +  {  (C1  -  C2 ) / b)x) sin  (nnx/b) dx 
b 

=  (2/b) j  (C  -  C.  )  sin  (nrrx/b)  dx 
0  u  l 

b 

+  (2/b)  /  (  (C,  -  C~) /b)xsin (nox/b) dx. 

0  1 

From  tables,  we  find 


/sinx  dx  =  -cosx  +  C 


/xsinx  dx  =  sinx  -  xcosx  +  C. 


Then  we  have 


C  =  -(2/mr)(C.,  -  C.  )  cos  (n~x/b) 
n  u  1 


2  iD  i° 

+  (2(C,  -  C_)/(mr)  )  (sin  (n~x/b)  l  -  (nTrx/b)  cos  (n:rx/b)  j  ) 

12  '0  '-0 
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r*y 

y 


C  =  (2/nrr)  (C,  (-1) n-C. )  +  (2Cn/nr)  (1-  (-l)n). 
n  2  l  u 


['ha  second  term  above  exists,  only  for  n  odd.  So  we  can 


define  it  as  4CQ/(2n  +  1),  n  =  0,1,2 .  Then  U  is  the 


form 


U  (x  ,  t )  =  (  2/tt)  y  {  (-l)n  (C,  -  C.  ) /n)  sin  (nTrx/b)  exp  (- {  (n-r)  "/ 

n=l  *  L 


b  j  Dt)  +  < 4 C _ / ?r )  l  (sin  (  (  (2n  -r  1)  /b)  -x)  / 


n=0 


(2n  +  1)  )  exp  (  -  (  (2n  +  1)  “fr‘*/b“)Dt) 


■.ere  tore 


C  { x  ,  t )  =  C1  +  ({C2  -  C1)/b: 


+  (2/tt)  £  (  (-l)nC2  -  C^)  /n)  sin  (nrx/b)  e 
n=l 


- ( (nr) “/b“) Dt 


+  (4Cn/TT)  l  (sin  (  ( 2n  +  l)/b)rx)/(2n  +  1) 


n=0 


exp(-(2n  +  1)  ^TT^/bx“ )  Dt) 


(90) 


The  quantity  of  gas  molecules  through  the  membrane  to  the 
downstream  of  the  chamber  per  second  is  then 


dQ/dt  =  DA(9C/9x)  where  Q  is  the  total  amount  of  diffusant 

JC  ~ "  u 


collected  in  the  downstream  reservoir,  D  is  the  diffusion 
coefficient  and  A  is  the  effective  cross  section  of  the 
~e.aorr.ne  normal  to  the  flow  direction  . 


The  total  amount  of  gas  collected  downstream  up  to 


time  t  i 


Q,  =  AD J  (3C/3X)  dt 
c  0  x  J 


Q_/ AD  =  /(  (C2  -  C1)/b)  +  (l/b)  l  ({-1  )nC2-C1) 
0  n=l 


-  (nrr/b)  ^Dt 


+  (4Cn/b)  J  e  /b  ;>'r  )Dt)dt 


n=0 


?  n 

=  {  (C9-Ci),/b}t  -  (2b/D:r)  T  (  (  ( -1 )  nC2  -  C. )  /n 

n=l 


(l-e“(n7r/b)  Dt)  -  (4Cnb/ba2)  £ 

n=() 


(l/(2n+l,2l(l-e-«l,2:’+1)2’2>/b2)Dt) 


(91 


Experimental  condizions  in  this  work  are  that  both  Cq  and 
are  zero,  i.e.  the  membrane  is  initially  evacuated  to 
zero  concentration  and  the  concentration  at  the  face 
through  which  diffusing  gas  molecule  emerges  is  maintained 
effectively  at  zero  concentration. 

Therefore  equation  (91)  becomes 

CO 

Qt/AD  =  ( C  -,  /  b )  t  -  ( 2 b / D t 2 )  T  (-l)nc9/n2 
w  "  n-1  Z 

+  ( 2b/D~2)  l  (  (-l)nC9/n2)e"D(nTT) 


I  (-1)  /n“  =  tt^/I  2 


Q J./.\C2b  =  Dt/b4*  +  (2/"“)  {7-/12) 


00  / 

2  o  ,.n,  2.  -D(nir/b)“t 
-  ( 2/tt  )  (  (-x)  /n  )  e 

n=l 

as  t  approaches  the  exponential  term  approaches  zero 
and  a  linear  approximation  results. 


Qt/A  a  (C2D/b) (t  -  b  /6D) 


The  interceot  of  cha  line  CL/A  vs  t  wioh  the  t  axis  is 
called  the  "time  lag." 


i.e.  t  =  b'VSD 


Therefore 


D  =  b  2  /  6  t 


If  Cn  and  2n  have  nonvanishing  values,  the  relation  is 


D  =  b 2  { C 2  +  2C1  -  3C0)/6t(C2  -  C^. 
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APPENDIX  C  MASS  UPTAKE  FOP.  UNIFORM 


SURFACE  CONCENTRATION 


The  boundary  conditions  are  below: 


C(x,t)  =  C?  ,  at  x  =  0,  all  t  >  0 


C(x,t)  =  C^,  at  x  =  b,  all  t  >  0 


C(x,t)  =  CQ,  at  0  <  x  <  b,  all  t  =  0 

The  only  difference  of  these  boundary  conditions  to 
equation  (83)  of  Appendix  3  is  C(x,t)  =  C^  at  x  =  b.  So, 
the  equation  (90)  of  Appendix  B  can  be  directly  applied  t 
here  by  changing  C2  to  C^. 

We  can  obtain 


C  ( x , t )  =  C,  +  (2/a)  /  ({(-i)nC1  -  C, ) /n) sin (nnx/b) 

1  n=l 

a-{(r.,)2/b2Dt)+  (4C  ,,,  | 

n=0 

(sin  (^~)  ax/ (2n+l)  )e  (  (2n+i)  71  /b  )Dt) 


The  second  term  in  above  equation  can  exist  only  n  is  odd 
nuiTioer.  The  equation  changes  to 


’mimw 


q-'. -  -'./.n''.'  ■;  v. ••vy-vV v;  ■  \ y. , 


C(x,t)  =  C,  +  ( 4 /  rr )  ( C „  -  C,  )  ]  (sin  (  (2n+l)  i'z  )  ax)  /  (2r.- .. 

1  n=0 

exp  (-  (  (2n+l)  2 7T 2 / b 2 )  Dt) 

by  changing  n  to  2n+l  in  the  second  terra  and  combining  with 
the  third  term.  After  rearrangement,  the  equation  becomes 


CO 

(C  -  -  C  )  =  1  -  (4/a)  y  ( 1  /  ( 2n  +  i)  ) 

u  x  u  n=0 

exp(-((2n  +  1)  ^TT^/b^)  Dt)  sin  (  (2n  +  l)/b)7TX  (93) 


We  define  it  as  the  amount  of  diffusant  has  diffused  into 
the  film  at  time  t,  and  Q  as  the  amount  of  diffusant  at 
infinite  time. 

b 

Because  Q  =  Jc(x,t)dx, 

0 

equation  (93)  can  be  simply  integrated  on  each  side  to  gat 


2  2  ^ 

Qot/  (C^-Cq)  =1-  (4/ti)  l  (1/ (2n+l)  )  exp  (- (  (2n+l)  /b  )  tt  Dt) 
b 

/sin  (  (2n  +  l) /b)  TTXdx  (94) 

0 

b 

Solve  ! sin (  ( 2n  +  l ) /b ) ax  dx, 

0 

we  can  get  2b/ (2n  +l)a 

this  back  to  equation  (94) ,  we  obtained 


( 8  /  ( 2n  +  l ) 


'xn  f-n  ■  ( 2 n + 1 )  ■“)  /b“)  t 


Put 


14  7 
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APPENDIX  D  CORRECTION  OF  EDGE  EFFECTS 

For  a  semi-infinite  solid  (h  >  0)  the  boundary  con¬ 
ditions  become 

C  =  C^,  0  <  x  <  “,  t^O 

C=C  ,  x  =  0 ,  t>0 
m 

So, 

(C  -  C.  )  /  (C  -  C.  )  =  1  -  erf  (x/2/D~t)  (95) 

X  III  1  .X 

and  the  weight  of  moisture  which  enters  or  leaves  through 
an  area  A  in  time  t  is 
t 

m  =  /-AD  (3C/3x)  dt  (95) 

0 

Equation  (95)  and  (96)  give 

m  =  2A(Cm  -  Ci)  /Dxt/ir  (97) 

In  the  early  stages  of  the  process  of  interaction  of  the 
different  sides  might  be  neglected  and  equation  (97)  can 
be  applied  to  each  side  independently.  Thus  the  total 
weight  of  entering  moisture  is 

m  —  4  ve  —  C  ■  )  (nl»D_  +  nh/D  +  nl»  D  )  v  t/  i 


JW.'M'IOT 


1 

v.'here  D  ,  D  ,  D  are  che  dif f usivities  in  the  x,  y,  z, 
x  y  z  >  i  >  > 

directions . 

Because  Ck  =  0 ,  =  oghni ,  the  mass  uptake  will  be 

M  =  (  (W  -  W,)/W.)  x  ICO  =  (m/W.)  x  100 
c.  a  d 

=  (4C  /ugh)  (/D~  +  (h/1)  v’D-  +  (h/r.) /D~ )  y  t/77 
m  ^  x  y  z 

C  / pc  is  the  maximum  moisture  content  as  in  the  ecru cion 
m  - 

shown  below. 


M  =  {  (V!  -  W,)/W.)  x  ICO 

m  m  d  d 

=  ( (m  /hnl)  /  (W./hnl)  )  x  100 

in  Cl 

-  cm/p9 


M  =  (4Min/h/:f }  {-.  Dx  +  (h/l)vD  +  (h/’n)^F)  t  c 

For  the  linear  part  in  the  beginning  of  a  diffusion  curve 
it  can  be  expressed  as 

M  =  (4M  /h)  yTT^v  D 
m 

then 


D 


(9 


APPENDIX  E  ESTIMATION  OF  PARAMETERS 


We  assume  10  centers  per  monomer  unit  in  PPBT,  then 

o  *  O 

X  equals  1.235  A.  For  p  we  choose  4.3  A  which  is  the 

O  O 

average  of  6.1  A  and  3.5  A  (two  transverse  distances  between 

46 

chains).  According  to  Small's  suggestion  that  the 
cohesive  energy  density  should  be  an  additive  property  for 
molecules,  similar  to  the  molar  refraction  or  the  parachor. 
5~of  PPBT  can  be  estimated  as  below: 


Molecular  weicht  for  PPBT  monomer: 


Density:  1.6 

1  phenylene  (p) 
1  phenylene  (m) 


C  =  N 


Ring,  5-membered 
Sul f ide 


6  =  d;:G/M  =  1.6  x  3146/267  =  18.9 
6^  =  355.4  cal/cm^  =  1487j/cm^ 


m  =  267/10  =  26.7  dalton 


Insert  above  values  into  equation 


=  : .  /(  3 . 88  (density ) .  ) 


3146 


r,rr?f 


ib  4ou  i  rule 


’■'3  obtain  £  = 

Because  a  molecule  ordinarily  vibrates  about  its 
equilibrium  position  at  a  frequency  of  lO1-  to  10l3 

:  2  -1 

vibrations  p^r  second,  we  assume  Vg  =  10  sec  .  Tri 
X  =£*p*/A2  =  1645  x  4.8  x  1(T8/  (1.235  x  10~8)2 
=  5.177  x  10^  Jcm/mole 

2  2 

D  =  vQ (RT)  /x  at  300°K 

-  1012 (8.314  X  300) 2 / (5.177  x  10"1)2 
=  2.32  x  1 0~D  cffl2/s 
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